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Although fundamental for astrophysics, the processes that produce 
massive stars are not well understood. Large distances, high extinction, 
and short timescales of critical evolutionary phases make observations of 
these processes challenging. Lacking good observational guidance, the- 
oretical models have remained controversial. This review offers a basic 
description of the collapse of a massive molecular core and a critical dis- 
cussion of the three competing concepts of massive star formation: 

• monolithic collapse in isolated cores 

• competitive accretion in a protocluster environment 

• stellar collisions and mergers in very dense systems 

We also review the observed outflows, multiplicity, and clustering proper- 
ties of massive stars, the upper initial mass function and the upper mass 
limit. We conclude that high-mass star formation is not merely a scaled- 
up version of low-mass star formation with higher accretion rates, but 
partly a mechanism of its own, primarily owing to the role of stellar mass 
and radiation pressure in controlling the dynamics. 
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1 INTRODUCTION 

1.1 Basic Issues 

Massive stars play a key role in the evolution of the Universe. They are the 
principal source of heavy elements and UV radiation. Through a combination 
of winds, massive outflows, expanding HII regions, and supernova explosions 

they provide an important source of mixing and turbulence in the interstellar 
medium (ISM) of galaxies. Turbulence in combination with differential rotation 
drives galactic dynamos. Galactic magnetic fields are generated, interacting 
with supernova shock fronts that accelerate cosmic rays. Cosmic rays, UV 
radiation, and dissipation of turbulence are the principal sources of heating in 
the ISM, whereas heavy elements found in dust, molecules, and in atomic/ionic 
form ultimately are responsible for its cooling. Massive stars thus profoundly 
affect the star- and planet-formation process (Bally, Moeckel & Throop 2005) 
as well as the physical, chemical, and morphological structure of galaxies (e.g., 
Kennicutt 1998, 2005). 

In spite of the dominant role that massive stars play in shaping galactic 
structure and evolution, our understanding of their formation and early evolu- 
tion is still sketchy. There are many reasons. High dust extinction makes it 
difficult to observe high-mass stars during critical early formation phases. They 
are rare. They evolve quickly and important evolutionary phases are short- 
lived. The theoretical problem is extremely complex. Finally, massive stars are 
seldomly (if at all) formed in isolation; the proximity of other high-mass stars 
compounds the complex influence of the forming star on its local environment 
via gravitational interactions, powerful outflows and winds, ionizing radiation, 
and supernovae. 

The low number statistics of young or forming high-mass stars is only par- 
tially offset by their higher luminosities, which allow us to study them at greater 
distances than their low-mass counterparts. However, insufficient spatial resolu- 
tion is an issue - an entire OB-star cluster is often contained in a single observing 
pixel (e.g., Henning & Stecklum 2002). 



6 



Zinnecker Sz YbrJce 



1.2 Definitions 

Star formation typically starts with a collapsing gas condensation (core) inside 
a larger subunit (clump) of a molecular cloud (cf. Williams, Blitz & McKee 
2000). A protostar forms that increases its mass by accretion (accumulation) 
of neighboring gas, while at the same time some mass loss occurs through a 
bipolar outflow and/or a coUimated jet. Let us define some of the terminology 
adopted here. 

One of the most misused terms in papers dealing with star formation is proto- 
star, which is considered the Holy Grail (Wynn- Williams 1982) of IR astronomy. 
Here, we reserve the term protostar or protostellar object for a gaseous object in 
hydrostatic equilibrium (gas pressure forces balance the gas' self-gravity), which 
has not yet begun hydrogen burning but which will, given time, burn hydro- 
gen. At the point hydrogen burning commences, we shall speak of a zero-age 
main-sequence (ZAMS) star; as long as hydrogen burning occurs in the center, 
we shall speak of a main-sequence star. Note that the size scale of a protostar 
is at most a few tens of solar radii. 

We use the terms massive star and high-mass star interchangeably to denote 
an OB star sufficiently massive to produce a type II supernova (M^/M© > 8 
for solax abundances). With these definitions in mind, the term high- mass pro- 
tostar denotes a > 8M0 hydrostatic object that has not yet begun hydrogen 
burning. As we shall see in the following, such objects exist only briefly during a 
transitory stage between "accreting intermediate- mass protostar" and "accret- 
ing high-mass star." Because it will be impossible to distinguish observationally 
when an accreting object begins burning hydrogen, we suggest that the terms 
massive protostar and high-mass protostar generally be avoided. 

In Table 1, we give a crude classification of massive stars in terms of loga- 
rithmic mass intervals and the corresponding main sequence spectral types. 

We reserve the terms very massive star (VMS) and supermassive star (SMS) 
for stars in the mass ranges of 100 <, M,/Mq ^ 1000 and 10* <, M^/Mq ^ 10^ , 
respectively, and introduce the term ultramassive star (UMS) for stars in the 
mass range of 10^ <^ M^/M© <^ 10*. SMSs are equilibrium configurations that are 
dominated by radiation pressure — baryons and electron-positron pairs provide 



Massive Star Formation 7 
Table 1: Main Sequence massive star definition (logarithmic mass ranges) 



Mass 


Designation 


Sp. type 


8-16Mq 


Early B-type massive stars 


B3V to BOV 


16-32Mq 


Late 0-type massive stars 


09V to 06V 


32-64Mq 


Early 0-type massive stars 


05V to 02V° 


64-128 Mq 


0/WR-typc massive stars 


WNL-H^ 



"02V main sequence stars have been identified by Walborn et al. (2002). 
''WNL-H: N-rich latc-type Wolf-Rayet (WR) stars, still on the 
Main Sequence (H-burning) - see Crowther (2007). 

only a minor contribution to the equation of state. At some point during their 

evolution SMSs collapse owing to a general relativistic gravitational instability. 
Whereas in the present epoch VMSs, UMSs, and SMSs are unlikely to be formed 
except under very special conditions, stars with masses in excess of 100 M© are 
expected during the first epoch of star formation (Bromm & Larson 2004; Abel, 
Bryan & Norman 2000). VMSs, UMSs, and SMSs arc not discussed in this 
review. A recent discussion of the formation and evolution of VMSs is given in 
Portcgics Zwart et al. (2006) and in Belkus, van Bever fc Vanbeveren (2007). 

What is accretion? The term accretion is used in a variety of senses. Mea- 
sured accretion rates often refer to the rate of mass inflow toward star-forming 
sites - not the rate at which a star or protostar gains mass. Originating from 
a 0.1 pc scale, this material cannot possibly fall into a sub-10~^ pc region with- 
out carrying signiflcant angular momentum. Instead, it either forms a disk or 
hits and is mixed with prior existing disk material (see Figure 1). Thus, we 
distinguish between the accretion of cloud core material onto a disk (Mn-acc) 
and the accretion onto a (proto-)star (Ms-acc)- 

Analogous to accretion, mass loss is used in a variety of senses. Measured 
mass loss rates from jets and outflows do not necessarily reflect the mass loss 
from an isolated young star. We thus distinguish between the mass loss from the 
(proto-)star via a wind (Ms-wind)) the mass loss launched from the accretion 
disk (MD-wind)) which never reaches the (proto-)star, and the material swept up 
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into the outflow from the surrounding molecular cloud (Mioad)- The measured 
outflow could have contributions from several stars and several disks. 

The interrelation between disk accretion and disk winds is a fascinating 
aspect of massive star formation (sec, e.g., the recent magneto- hydrodynamic 
models of Bancrjcc & Pudritz 2007) and is at the focus of the frequently 
asked question: Is high-mass star formation a scaled-up version of low-mass 
star formation? The answer to this question will be better defined by the time 
we reach the end of this review. 

1.3 Recommended Reading 

The study of the origin of massive stars is a relatively new field of astrophysical 
research. There is no comprehensive monograph on the subject, but there are 
several conference proceedings over the past few years dedicated to the topic, 
of which we recommend the following; 

• Massive Stars: Their Lives in the interstellar Medium (Cassinelli & Church- 
well 1993) 

• Hot Star Workshop III: The Earliest Stages of Massive Star Birth (Crowther 
2002) 

• Massive Star Birth: A Crossroads of Astrophysics (Cesaroni et al. 2005a) 

The reviews on Environment and Formation of Massive Stars (Garay & 
Lizano 1999), Control of Star Formation by Supersonic Turbulence (Mac Low 
& Klesscn 2004), The Formation of the First Stars in the Universe (Glover 
2005), The Birth of Massive Stars and Star Clusters (Tan 2005), High Mass 
Star Formation by Gravitational Collapse of Massive Cores (Krumholz 2007), 
and The Critical Role of Disks in the Formation of High-Mass Stars (Cesaroni 
et al. 2007) are also recommended. Among the reviews in the Proceedings of 
Protostars and Planets V, we particularly recommend The Formation of Massive 
Stars (Beuther et al. 2007). 

There are also a few related Annual Reviews articles: 

• Compact HU Regions and OB Star Formation (Habing & Israel 1979) 
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• The Search for Infrared Protostars (Wynn- Williams 1982) 

• The Dynamic Evolution of HII Regions - Recent Theoretical Developments 
(Yorke 1986) 

• The Orion Molecular Cloud and Star-Forming Region (Genzel & Stutzki 
1989) 

• Physical Conditions in Regions of Star Formation (Evans 1999) 

• Ultra- Compact HII Regions and Massive Star Formation (Churchwell 2002) 

• Massive Stars in the Local Group: Implications for Stellar Evolution and 
Star Formation (Massey 2003) 

• The First Stars (Bromm & Larson 2004) 

Except for the last, none of these earlier reviews had its focus on the for- 
mation aspect of massive stars but rather provided a descriptive observational 
summary of the properties of young OB stars and their HII regions. The present 
Annual Reviews article will be accompanied in the same volume by an article 
on 'Physical properties of Wolf-Rayet stars' (Crowther 2007) and an article on 
'Theory of Star Formation' (McKee & Ostriker 2007) which mostly addresses 
low-mass star formation, but also includes an important section on high-mass 
star formation. 

1.4 The Focus of This Review 

The major questions we wish to address in this review are: 

1. What is the sequence of observable states leading from molecular clouds 
to young high-mass stars? 

2. What arc the initial conditions of massive star formation (gas densities, 
temperatures, clump masses, etc.) and how do they come about? 

3. Do massive stars always form in dense stellar clusters or can they form in 
isolation? What special conditions are necessary to allow coalescence, i.e., 
mergers of stars? 
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4. Which clues to the origin can be gleaned from multiplicity observations? 
How do we explain the very tight massive spectroscopic binaries and OB 
runaway stars? 

5. How does the forming massive star influence its immediate surroundings, 
possibly limiting its final mass and/or the final mass of its neighbors? 

6. How do young massive stars influence their global environment, either 
by inhibiting or by triggering further star formation? How do we get a 
starburst? 

To tackle these questions, we first discuss some key observations related to 
massive star formation. 
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2 MASSIVE STAR FORMATION: 
KEY OBSERVATIONS 

2.1 Observable Stages 

The optically visible main-sequence life of OB-type stars is preceded by an 
embedded phase that lasts about 15% of their lifetime (Churchwell 2002). As 
summarized by Menten, Pillai & Wyrowski (2005) and van der Tak & Menten 
(2005), observations at mid-IR through radio wavelengths have shown that this 
embedded phase can be subdivided into several groups of objects: 

• IR dark clouds (Perault et al. 1996, ISOCAM; Egan et al. 1998, MSX; 
Benjamin et al. 2003, Spitzer). Their internal density maxima and tem- 
perature minima likely represent the initial conditions of high-mass star 
formation; a compilation of several dozens of such high-mass starless cores 
has been given by Sridharan, Williams & Fuller (2005). Some of these 
cores probably contain low-mass and intermediate-mass accreting proto- 
stars, which are faint and hard to detect. Protostellar outflow activity has 
been detected in one of them (Beuther, Sridharan & Saito 2005). 

• Hot molecular cores (Kurtz et al. 2000, Cesaroni 2005). These have large 
masses of warm and dense gas, and large abundances of complex organic 
molecules evaporated off dust grains; they are signposted by methanol 
maser emission (Menten 1991, Walsh et al. 1998, Hill et al. 2005); 
ground-based dctectability on the Wicn part of the spectral energy distri- 
bution with sufficient spatial resolution is difficult (Stccklum et al. 2002), 
but comes into reach with dedicated 8-m class telescope observations (De 
Buizer & Minier 2005, Linz et al. 2005). 

• Hypercompact and ultracompact HII regions (Kurtz 2005, Hoare et al. 
2007). In these regions, small but growing pockets of ionized gas have de- 
veloped that stay confined to the stellar vicinity. Whereas hypercompact 
HII regions probably represent individual photoevaporating disks (Keto 
2007; see also the example in Nielbock et al. 2007), ultracompact HII 
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regions probably represent disk-less stars photoionizing their own cocoons 
and massive envelopes. 

• Compact and classical HII regions (Mezger et al. 1967, Yorke 1986). Their 
gas is ionized globally, often by several ionizing sources. It expands hydro- 
dynamically as a whole and disrupts the parent molecular cloud, revealing 
both the embedded high-mass and lower mass stellar population for optical 
and near-IR observations (Carpenter et al. 1993; Zinnecker, McCaughrean 
& Wilking 1993). 

2.2 Initial Conditions 

Massive star formation occurs inside dense, compact clumps in giant molec- 
ular clouds (H2 column densities are 10^^-10^^cm~^). Smaller mass clumps 
with lower peak H2 column densities do not form massive stars. Several types 
of molecular cloud surveys, predominantly near HII regions, have been carried 
out: CS-molecule surveys for dense molecular gas, 1.2-mm dust continuum sur- 
veys for massive cold dust (and hence gas) condensations, as well as OH, H2O, 
and methanol maser emission surveys for shock-excited compact regions as sign- 
posts for massive star formation. (Note that methanol maser and OH maser 
emission is exclusively associated with high-mass star formation, whereas H2O 
mascrs may also bo found in low-mass star-forming regions. This is because 
methanol and OH mascrs arc radiativcly pumped and need an intense far-IR 
source in their vicinity; H2O mascrs, in contrast, are coUisionally pumped in 
gas shocked by outflows.) These gas and dust surveys have revealed dense cold 
clumps (molecular hydrogen density = 10^cm~^, gas temperature T = 10- 
20 K, diameter -^O.Spc) with gas masses ranging from a few hundred to a few 
thousand solar masses (Plume et al. 1997, Shirley et al. 2003, Garay et al. 
2004, Motte et al. 2005, Evans 2005). The methanol maser surveys at 6.7 GHz 
point to hot molecular cores with internal heat sources and outflows, as well as 
protoclusters (Burton et al. 2005, De Buizer 2003, Minier et al. 2005). 

As mentioned before, large-scale observations with the ISO, MSX, and Spitzer 
satellites have revealed a new class of clouds, the so-called IR dark clouds or 
IRDCs, and Simon et al. (2006) have identified more than 10,000 such IRDCs 
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from the MSX data base. Many of these appear to be located in the 4-5 kpc 

Galactic molecular ring (Jackson, private communication; see also Bronfman 
et al. 2000). The IRDCs arc dense clouds seen in absorption against mid-IR 
background emission. They arc mostly filamentary structures that contain con- 
densations of cold massive cores where massive stars or even star clusters seem 
to form (Rathbornc, Jackson & Simon 2006). Recent mid-IR and millimeter- 
continuum observations show different evolutionary stages of massive star for- 
mation in adjacent cores: dense millimeter-continuum sources with and without 
mid-IR emission (Garay et al. 2004, their figure 4). 

The origin of these structures appears to derive from supersonic turbulence 
in giant molecular clouds, that is, shock compression from convergent turbu- 
lent gas streams. Depending on the direction of the compression with respect 
to the direction of the magnetic field lines, the magnetic field will be boosted 
through flux freezing and hence the resulting clump will be stabilized by mag- 
netic forces against gravitational collapse (subcritical compression). If not, the 
compressed clump is quickly set up for collapse (supercritical compression), in 
fact so quickly that the set-up time is shorter than the free-fall time. Mag- 
netically stabilized clumps take much longer for collapse to begin, and their 
internal turbulent structure may make the clump prone to subfragmentation. 
If this is true, only supercritical compression leads to massive star formation 
(cf. Shu, Adams & Lizano 1987). An interesting speculation would be that 
clouds without sufficiently strong magnetic fields can form lots of massive stars 
quasi-simultaneously, giving rise to gigantic starbursts. 

2.3 Endproducts 
2.3.1 OB Clusters 

The endproducts of massive star formation are either dense gravitationally 
bound OB star clusters or loose unbound OB associations (Lada & Lada 2003, 
Briceno et al. 2007). Classical examples of OB star clusters include the Orion 
Nebula Cluster (ONC), the dense compact cluster associated with the giant 
galactic HII region NGC 3603, and the R136 cluster in the 30 Dor region in the 
Large Magellanic Cloud (LMC; see Figure 2). These clusters roughly define 
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a richness sequence in powers of 10: they contain 1, 10 (21), and ~ 100 mas- 
sive 0-type stars per cluster, with estimated total cluster masses of 10^, 10^, 
and 10^ M0, respectively (e.g., ONC: Hillenbrand 1997, Hillenbrand & Hart- 
mann 1998; NGC 3603: Moffat, Drissen & Shara 1994, Drisscn ct al. 1995; 
R136: Parker & Garmany 1993, Masscy & Hunter 1998). Although R136 can 
probably be considered a small young globular cluster (M. Andersen ct al., sub- 
mitted), there are more massive young clusters in the nearby universe, such as 
the very massive embedded super star cluster in the center of the NGC 5253 
dwarf galaxy (Turner et al. 2003). This cluster has an ionizing flux equivalent 
to the presence of 4000-6000 07V stars that comes from a very compact region 
about 1 pc in size, measured with the VLA. How can the formation of so many 
massive stars in such a small volume be possible? This question is at the heart 
of understanding the origin of globular clusters. 

2.3.2 OB Associations 

The classical examples of OB associations include the nearby Scorpius 0B2 and 
Orion OBI associations (Blaauw 1964, 1991). Another fine example of what 
may ultimately become an OB association is the Carina star formation complex 
at 2.3 kpc (Smith & Brooks 2007). One of the best studied extragalactic OB 
associations is NGC 604 in M33 (Mafz-Apellaniz, Perez & Mas-Hesse 2004). 
In all these cases, the OB stars are spread over the whole face of the parent 
giant molecular cloud and are not densely packed at all, with distances between 
massive stars ranging from 1 to 10 pc. This then appears to be a completely 
different mode of massive star formation, although it must be noted that OB 
associations often contain dense clusters, too (e.g., the Carina complex harbors 
the well-known Trumpler 14 and 16 clusters). The question is whether OB as- 
sociations arc superpositions of expanded young clusters (e.g., Kroupa, Aarseth 
& Hurley 2001, Bastian & Goodwin 2006). 

2.3.3 Field OB Stars 

Do massive stars always occur in young star clusters or OB associations, or can 
massive stars also be found outside these regions, i.e., in the field? The answer is 
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they can. It has long been reaUzed that there exists a class of massive stars, the 
so-called runaway OB stars (Blaauw 1961, Poveda, Ruiz & Allen 1967, Gies & 
Bolton 1986) that are ejected from their birthplaces - clusters and associations 
with velocities in excess of 40kms^^. About 10-25% of all O stars and about 
2 % of all B stars belong to this class of massive field stars. The question whether 
- these runaway stars aside other massive stars occur in the field (implying 
that they would be formed in isolation) has been studied by de Wit et al. (2004, 
2005), following Mason et al. (1998). These authors established that 43 among 
the 227 O stars brighter than eighth V-magnitude are in the field. Of these, 
about half can be traced back to a cluster or association origin, but about 10 - 20 
stars (i.e., 5 - 10 %) could not be assigned to any group of origin and therefore 
might be true field stars, born outside clusters and associations, an issue already 
raised in a pioneering paper by Roberts (1957). A case in point is HD93521, 
a high-latitude 09.5V star, more than 1 kpc above the Galactic plane, which 
must have formed locally in the halo (Irvine 1989)! These examples suggest that 
the question of the birthplaces of massive stars is not yet completely settled. 
The forthcoming Spitzer 8- micron imaging survey of the Magellanic Clouds (see 
Mcixncr et al. 2006) can shed new light upon the cjucstion, and in particular can 
pinpoint isolated massive stars in the LMC/SMC, should these objects indeed 
exist. 

2.4 Clues from multiplicity 

The multiplicity of massive stars is believed to be higher than that of young low- 
mass premain-sequcnce stars (Preibisch, Weigelt & Zinnecker 2001, Duchene 
et al. 2001). This means massive stars have more physical companions than 
low-mass stars on average. For reference, the multiplicity or, more precisely, 
the companion star fraction (csf) of a stellar population has been defined by 
Reipurth & Zinnecker (1993) to be 

csf = (B -h 2T + 3Q + ...)/(S + B -h T + Q + ...), 

where S is the number of single, B the number of binary, T the number of triple, 
and Q the number of quadruple systems, etc. (i.e., triple systems contribute 
two companions, quadruple systems three companions, etc.). For example, the 
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multiplicity of the four OB stars in the Trapezium Cluster in Orion is as follows: 

the most massive star &^ C is double, the next massive star A is triple (a 
hierarchical system with a close spectroscopic binary and a wider companion), 

B is at least quadruple (kind of a Trapezium system within the Trapezium 
Cluster), and D is apparently single. [Kraus et al. (2007) find some indica- 
tions that 01 D appears extended in their speckle images. There is a fifth star 
in the Trapezium Cluster, E, which has recently been found to be a double- 
lined intermediate-mass spetroscopic binary (Herbig & Griffin 2006, Costero et 
al. 2006).] Putting the number of companions of A, B, C, D into the above 
formula, we get csf=1.5 at face value (probably a lower limit). This should 
be compared to the multiplicity of the low-mass stellar members in the Orion 
Nebula cluster, which is csf=0.5 (Padgett, Strom & Ghez 1997; B. Reipurth et 
al., submitted), i.e., significantly lower. 

The multiplicity statistics of other OB clusters, including clusters rich in O 
stars (N > 5) and poor in O stars (N < 3) has been studied by Mermilliod & 
Garcia (2001) and Garcia & Mermilliod (2001), with interesting results: The 
spectroscopic binary frequency in 0-star rich clusters can vary enormously in 
different clusters, from 15% to 80%, with no apparent correlation. If anything, 
there is an anticorrelation of the binary frequency and the cluster density, but 
this needs to be reinvestigated and confirmed. The above statistics often rely on 
relatively poor data; and the sampling is not complete. Recent work on higher 
quality data, but on a more limited number of clusters, tends to obtain a lower 
binary fraction in the range of 20 % to 60 % (Sana, Rauw & Gosset 2005). The 
most dramatic example is the IC 1805 cluster where the binary frequency went 
down from 80 % to 20 % based on better data (De Becker et al. 2006). 

In the 0-star poor clusters almost all O stars are spectroscopic binaries, 
often double-lined and even eclipsing. These massive binaries are usually mem- 
bers of hierarchical triple or quadruple systems, or of trapezia, and are often 
located at the cluster center. The exciting star of the Orion Nebula Cluster, 

Ori C belongs in this category, although Ori C is not a massive close 
spectroscopic binary but a very eccentric visual binary, with masses of 34 
(05.5V) and 15 (09.5V) and an orbital period of about 11 yr (Kraus et 
al. 2007). The orbital periods of the spectroscopic binaries in the 0-star rich 
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clusters are concentrated in the range of 4-5 days, whereas in the 0-star poor 
clusters there is a pile-up of orbital periods around 3±1 days. In NGC 6231, 
for example, according to Sana (private communication), 10 out of the 16 O 
stars are double-lined spectroscopic binaries: 6 with periods under 10 days (4 
below 5 days); 2 with periods between 3 and 9 months, and 2 with periods of 
the order of a year or greater (in addition, 1 star is probably a triple-lined spec- 
troscopic binary). The luminosity ratios are all in the range 1-10, otherwise 
one would not detect them as double-lined spectroscopic binaries. This implies 
that the secondaries are probably early B stars, and the primary-to-secondary 
mass ratios must be about 3 (at maximum) or lower. It is also worth mention- 
ing that there is no very short-period highly eccentric 0-1-0 binary known at 
this time (Sana, private communication). These fascinating and surprising facts 
challenge our views of massive star formation and provide clues to their origin, 
clues too complex to fully decipher yet but hinting at gravitational dynamics 
playing a role - beyond mere disk or filament fragmentation (see Section 5, 
where the formation of binary and multiple systems is extensively discussed; see 
also Zinnecker 2003). 

Of course, it is equally important to study the multiplicity of massive stars in 
OB associations where the stellar density is much lower than in OB clusters, and 
dynamical interactions between the forming massive stars should be less of an 
issue. For example, the Orion OBI association contains ^ 70 massive stars in its 
three subgroups la,b,c (subgroup Id is the Orion Nebula Cluster). Of these 70 
OB stars, 20 % are spectroscopic binaries with periods less than 10 days (Morrell 
& Levato 1991). The three subgroups show variations in their spectroscopic 
binary fractions: subgroup la is average, subgroup lb is a factor 1.5 above the 
average, and subgroup Ic is a factor 1.5 below the average. The situation in the 
nearby Scorpius-Centaurus 0B2 association is as follows: Among the 48 early 
B stars (B05 to B3V) there are 25 binaries, and 20 of them are spectroscopic 
binaries with known periods in the range of 0.9-34.2 days, with a median of 
5.7 days (Brown 2001). 

In conclusion, it seems that the spectroscopic binary fraction among massive 
stars in OB associations is surprisingly similar to that in OB clusters rich in O 
stars (about 40% on average), and the fraction of very close spectroscopic and 
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eclipsing massive binaries in OB associations with orbital periods below 5- 

10 days (about 20%, judging from the Orion and Scorpius-Centaurus regions) 
is a factor of two lower than in OB chisters. Thus OB clusters appear to 
contain more of a population of very tight (hard) binaries, possibly an effect 
owing to dynamical encounters after birth; this is an effect that is absent in OB 
associations. 

The future of spectroscopic massive binary research hes in the near-IR and 
in multiepoch radial velocity surveys of embedded massive stars. First results 
(and successes) have been reported by Apai et al. (2007), indicating that massive 
close binaries indeed form at a very early stage. 

2.5 Upper Initial Mass Function and Upper Mass Limit 

The mass distribution function of massive stars at birth (the so-called Initial 
Mass Function or IMF for short) is a complicated matter, because (a) massive 
stars quickly lose some of their initial mass through stellar winds, (6) many of 
these massive stars are unresolved binaries, and (c) massive stars tend to be born 
in the centers of OB clusters or tend to sink preferentially toward the cluster 
center, leaving behind their lower mass siblings that live in the cluster outskirts. 
This introduces a bias into the mass distribution, flattening a power-law slope. 
The upper IMF from about 10 to 100 Mq is usually found to be a universal 
power law, with logarithmic slope -1.35, first found by Salpeter (1955) for a 
range of masses below 10 M©. We refer here to the early observational work 
of Garmany, Conti & Chiosi (1982) and the summary of Massey (1998). The 
implication of a Salpeter slope or other similar slopes of the IMF for the number 
of stars born in different mass intervals (for convenience spaced by a factor of 
two) can be seen in Table 2, which has been normalized to contain exactly 
one object in the highest mass interval. This is instructive, because it shows 
dramatically how rare the 0/WR-type massive stars (interval 64-128M0) are 
compared with the early B-type massive stars (8 - 16 Mq), or with the solar-type 
low-mass stars (I-2M0); see Zinnecker (1996). 

The IMF of massive stars in the aforementioned OB clusters (Orion Nebula 
Cluster, NGC 3603, and R136 in 30 Dor) is discussed by Pflamm-Altenburg & 
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Table 2: Initial Mass Function (clN/dlogM ~ M"'') examples 



Mass range 




Logarillimic slojjo 






x=l 


x=1.35 


x=1.7 


0.5-1Mq 


128 


700 


3822 


1-2M0 


64 


275 


1176 


2-4Mq 


32 


108 


362 


4-8Mq 


16 


42 


111 


8-I6M0 


8 


16.6 


34.3 


16-32Mq 


4 


6.5 


10.6 


32~64M0 


2 


2.55 


3.25 


64-128Mq 


1 


1 


1 



Kroupa (2006), Stolte et al. (2006), and Massey & Hunter (1998), respectively. 
[Other studies of the Orion Nebula Cluster include Zinnecker, McCaughrean & 
Wilking (1993); Hillenbrand & Hartmann (1998). The stellar content in the 
biggest galactic HH region NGC 3603 was also investigated by Drissen et al. 
(1995); Hofmann, Seggewiss & Weigelt (1995); Eisenhauer et al. (1998); Brandl 
et al. (1999): Moffat et al. (2004); and Sung & Bessell (2004). An important 
investigation of the IMF of R136 is that of Sirianni et al. (2000).] Whereas 
Pflamm-Altenburg & Kroupa (2006) find a deficit of high-mass stars in the 
Orion cluster, Stolte et al. (2006) derive an excess of massive stars in the NGC 
3603 cluster core, reflected by a power-law slope of -0.9 (probably due to mass 
segregation). It is only in R136 that the power-law slope is almost exactly the 
same as the Salpeter value. It is worth noting here that the slope of the mass 
function of high-mass stars in the range of 8 - 40 Mq in the wider field in the 30 
Dor region is apparently the same as in the R136/NGC 2070 cluster (Selman 
& Melnick 2005). The latter authors do not find the much steeper slope of the 
IMF (for the range of 25 Mq to 120 Mq) derived by Massey (2002) for the global 
OB field population in the LMC. They suspect that selective incompleteness at 
V = 12 owing to detector saturation and Be star contamination lies at the origin 
of this discrepancy. 
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The question of the IMF of massive stars in OB associations was discussed 

long ago by Garmany, Conti & Masscy (1980) and later by Masscy, Johnson & 
Dcgioia-Eastwood (1995). They concluded that the IMF is normal, i.e., con- 
sistent with a Salpeter power-law. A fine discussion of upper IMF slopes in 
various young clusters and associations including all the caveats and selection 
effects was given by Scalo (1998). He noted that individual realizations of IMF 
slopes can vary, but the average slope is indeed close to the Salpeter value -1.35. 

The question of whether there is a physical (rather than statistical) upper 
mass end to the IMF is of great interest to anyone interested in population 
synthesis, galactic evolution, and cosmology. For example, if the stellar upper 
mass limit were 120 Mq, pair instability supernovae requiring stellar masses 
at the time of explosion between 140 and 260 M© (Heger et al. 2003) could 
not happen. Larson (1982) originally asked the observational question whether 
there was a correlation between the mass of a molecular cloud and the maximum 
mass of a star that could form in it. His result was that indeed the maximum 
stellar mass scaled with the mass of the parent cloud, roughly with its square 
root. It takes a giant molecular cloud of 10^ M© to form a 50 Mq star; a cloud 
of 10^ Mq can only spawn a maximum stellar mass of 8M0. The implication 
is that massive stars form in clouds of mass between 10'^ and 10''' Mq or more, 
probably because only these have sufficiently massive substructure (clumps). 

Weidner & Kroupa (2004) and Figer (2005) discussed the upper limit to 
the masses of stars, based on observations of the R136 cluster in the LMC 
and the Arches cluster near the Galactic Center, respectively. They pointed 
out that these clusters are so massive that given a Salpeter IMF one would 
expect to find stars as massive as 750 M© and 500 M©, respectively, whereas the 
most massive stars seen do not exceed 140 M© and 130 M©, respectively. This 
suggests a firm upper mass limit of 150 M©. Otherwise, a sharp down-turn of the 
IMF near 150 M© would be required (see the extensive discussion in Elmegreen 
2000). Or the very massive stars have already exploded/imploded during the 
dust-obscured, hidden, early evolutionary stages - an unlikely scenario. Oey & 
Clarke (2005) also gave a statistical confirmation of a stellar upper mass limit 
around 120- 200 Mq, if the IMF is Salpeter-like. Koen (2006) further analyzed 
the upper IMF in the R136 cluster with two different statistical techniques and 
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suggested an upper mass limit of 140- 160 Mq. Thus all four studies agree on 
the existence of a physical upper limit in the stellar mass distribution. 

2.6 Feedback and Triggering 

This topic deserves its own review. The question we ask here is the follow- 
ing: What does the energy and momentum input of massive stars in terms of 
expanding HII regions, stellar winds, or supernova shock waves do to the par- 
ent clouds? Is the cloud primarily disrupted or is new star formation triggered? 
Which of the above agents (HII regions, stellar winds, or supernova shock waves) 
provides the best trigger for new OB star formation and for new low-mass star 
formation? We arc only beginning to answer these questions. 

A key observation in this context is the fact that the high-mass and low-mass 
stellar populations in the subgroups of OB associations appear to be coeval, 
i.e., the nuclear age of the massive stars is the same as that of the lower mass 
premain-sequence objects (Preibisch & Zinnecker 1999, 2007; Bricefio et al. 
2007). This would appear to require a fast, coherent trigger, such as a supernova 
shock wave. Indeed, in the Scorpius-Centaurus association there is evidence 
that the shock wave of a supernova in one OB subgroup triggered the formation 
of another subgroup (de Geus 1992). However, there is also other evidence 
that radiation from massive stars [by a process called radiative implosion (e.g., 
Kessel-Deynet & Burkert 2003)] can only trigger the formation of low- and 
intermediate- mass objects (Lee & Chen 2007 and references therein). 

The classic theory of triggered massive star formation is the one by Elmegreen 
& Lada (1977). In their theory, an ionization shock front provides the pressure 
on an adjacent layer of molecular gas to compress it and heat it, thus stim- 
ulating gravitational instability of massive gas layers, hence the formation of 
massive stars - and perhaps only massive stars. Low-mass star formation may 
not be triggered with this mechanism. If true, this would produce an anomalous 
IMF with only massive stars (known as bimodal star formation, cf. Glisten & 
Mezger 1982). It is worth noting that such characteristics - small groups of 
massive stars only - may be needed for the decay of small N-body systems that 
give rise to the dynamical ejection of runaway OB stars (Clarke & Pringle 1992). 
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3 MASSIVE STAR FORMATION: 
BASIC THEORY 

3.1 Sequence of Events 

Starting from a pre-existing giant molecular cloud, the sequence of events is 
likely as follows: 

1. Formation of cold dense molecular cores or filaments, induced by gravo- 
turbulcnt cloud fragmentation (Mac Low & Klcsscn 2004). This means 
that supersonic turbulence rapidly produces localized compressed pockets 
of gas, some of which remain gravitationally bound and provide the initial 
conditions for collapse (Padoan & Nordlund 2002; Klessen et al. 2005). 
A characteristic density of about 10^ cm"^ and temperature of 10-15K 
results from the equation of state of dusty molecular gas of solar abundance 
(Jappsen et al. 2005, Larson 2005). 

2. Nonhomologous gravitational collapse of portions of the cores into opti- 
cally thick, pressure-supported protostellar embryos with initial masses of 
the order of 10^"^ (Larson 1969, Bate 2000). The term nonhomologous 
collapse refers to the fact that the relative distribution of material changes, 
as opposed to a homologous or self-similar collapse. (This is different from 
Shu's (1977) inside-out collapse of a self-similar isothermal sphere). 

3. Accretion of material onto protostellar objects as they evolve toward the 
main sequence. For low-mass objects the accretion stops well before hy- 
drogen burning commences. These premain-sequence objects of fixed mass 
then slowly and quasi- hydrostatically contract to the main sequence (Palla 
& Stahler 1993, Baraffe et al. 2002). However, high-mass objects even- 
tually start burning hydrogen and develop radiation-driven winds as they 
continue to accrete and evolve up the main sequence to hotter and more 
luminous states (Kudritzki 2002). 

4. Disruption of the birth cloud, as the first high-mass stars strongly influ- 
ence their environment by their winds, outflows, and UV radiation, and 
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eventually become supernovae. The most massive stars go supernova af- 
ter ~ 3 Myr. When the remnant molecular cloud has been dissipated, the 
result is mostly a chistcr of OB stars or an OB association, with an asso- 
ciated cospatial population of lower mass stars (Zinnecker, McCaughrean 
& Wilking 1993). Often, several evolutionary stages of star formation can 
be found side by side; the 30 Dor region in the LMC is a good example 
(Walborn et al. 1999). 

We shall denote these four phases as the compression, collapse, accretion, 
and disruption phases of high-mass star formation, respectively. All phases can 
occur simultaneously and sidc-by-side in a molecular cloud. Below, we give a 
fairly detailed yet still greatly simplified discussion of these phases, inchiding 
protostellar luminosity and ionization evolution as a function of growing stellar 
mass. 

3.2 The Compression Phase 

This first step toward high-mass star formation is either a starless core 100 Mq) 
or a starless clump 1000©) of molecular gas in a giant molecular cloud. Mc- 
Kee & Tan (2003) envisage that these cores are molecular condensations in a 
turbulence-supported quasi-equilibrium that ultimately form single or gravita- 
tionally bound multiple massive protostars. As these authors argue, turbulent 
and pressurized clouds permit sufficient material to be available in the cores of 
giant molecular clouds for high-mass star formation. Mechanical energy must be 
continuously injected into the clump in order to maintain this quasi-equilibrium 
between turbulence and gravity. The assumption is that this energy is either 
injected from within the cores from the kinetic energy of outflows and accretion 
shocks, or it comes from the outside and cascades down to smaller scale sizes. It 
is furthermore assumed that (small-scale) turbulence acts as an isotropic pres- 
sure. 

In an alternate scenario by Bonnell et al. (1997, 2001a), the compression 
phase is more of a transient phase due to the random motions in the self- 
gravitating cloud. Smoothed particle hydrodynamics (SPH, see the sidebar 
Smoothed Particle Hydrodynamics with Sink Particles) simulations of this phase 
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show sheets, filaments, and cores forming (see Figure 3). Some cores collapse 
and fragment, marking the beginning of the collapse phase for these objects, 
while in other parts of the cloud compression is still occurring. 

The philosophical difference between the approach of McKee & Tan (2003) 
and that of Bonnell et al. (1997, 2001a) can be summarized as "monolithic 
collapse" versus "competitive accretion." For the former, the mass necessary 
for massive star formation is intimately associated with the final product. If 
there are bulk motions of the embryo star, its protostellar core participates in 
those motions. The only competition for the infalling material is between close 
members of a multiple system. For competitive accretion, the material that 
makes up a particular star can come from various parts of the parent cloud. 
Protostars move relative to the molecular gas; the only gas that is intimately 
associated with any particular protostar is in its circumstellar disk and envelope. 
Because multiple protostars are often formed together at the same time, each 
protostar competes for the available molecular material (see Section 4.2). 

Crutcher's (1999, 2005) summary of the available Zeeman measurements of 
magnetic field strengths in molecular clouds suggests that magnetic fields likely 
play an important role in molecular cloud dynamics. Vazquez- Semadeni et al. 
(2005) studied the evolution of clumps and cores formed as turbulent density 
fluctuations in nearly isothermal molecular clouds, considering both the mag- 
netic and nonmagnetic cases with driven turbulence. In the nonmagnetic case 
the authors find that the cores arc unlikely to reach a hydrostatic state neces- 
sary for monolithic collapse - if the molecular clouds have an effective polytropic 
exponent less than 4/3. In this case, cores are transient, either proceeding di- 
rectly to collapse or re-expanding on a dynamical timescale. 

The magnetically subcritical clouds simulated by Vazquez-Semadeni et al. 
(2005) do not produce magnetostatic clumps, but rather a few marginally bound 
clumps that are subsequently dispersed. Ambipolar diffusion - had it been 
included in the simulations - could have increased the clumps' likelihood to 
become bound and subsequently to collapse. For clouds with weaker magnetic 
fields a few cores form and collapse on a timescale slightly larger than the 
cloud's free-fall timescale. In their most supercritical simulation, fewer clumps 
and cores form than in their nonmagnetic counterpart, and these cores reexpand 
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because they are not Jeans unstable. The authors thus conclude that not all 
cores observed in molecular clouds will necessarily form stars and that magnetic 
fields may help reduce the star-formation efficiency by reducing core formation 
rather than by delaying or inhibiting the collapse of individual cores. 
Smoothed Particle Hydrodynamics with Sink Particles 
Many simulations of star formation in turbulent interstellar clouds are conducted 
using a Lagrangian method to solve the equations of hydrodynamics. This is 
called smoothed particle hydrodynamics (SPH). It is a 3D method that assumes 
no symmetries and uses no imposed grids. Monaghan (2005) gives a recent 
review of the method and its applications (see also Monaghan 1992). With 
SPH, the gas is represented by a number of sampling points ("particles"), each 
associated with a certain gas mass. The mass of each particle is distributed in 
space according to a spread function so that density and other hydrodynamic 
quantities are continuous in space. The particles are allowed to move in the 
computational domain under the laws of self-gravitating hydrodynamics (which 
are evaluated at the particle positions). The system is evolved in time using a 
standard integration routine. 

The timestep with which the simulations evolve is chosen such that all 
changes to the system's state are examined at the appropriate speed. For ex- 
ample, during the collapse phase, when infalling particles approach each other, 
the timestep gets progressively shorter. This may potentially bring the sim- 
ulations to an early halt as the timestep converges to zero. This problem is 
overcome with the introduction of star (sink) particles of finite radii that each 
replace groups of dense neighboring gas particles above a certain density thresh- 
old (see Bate, Bonnell & Price 1995). The sink particles are then evolved as 
collisionless matter that can accrete from the lower density gas whenever gas 
particles approach within the sink radius. The final state of the system in such 
simulations involves groups of sink particles in spatial distributions resembling 
observed young clusters. However, all information on scales smaller than the 
sink radius is lost; such sink particles cannot later become sources of outflows, 
unless specifically instructed to do so. Note that sink particles are also usefully 
employed in Eulerian collapse calculations, as shown by Krumholz, McKee & 
Klein (2004). 
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In order for fragmentation to be realistically modelled, the numerical reso- 
lution of the SPH simulations must be such that they obey the Jeans condition. 
This states that "the local Jeans mass must always remain resolved," i.e., it must 
contain more than a certain number of particles (Truelove et al. 1997, Bate & 
Burkert 1997, Kitsionas & Whitworth 2002, Martel, Evans & Shapiro 2006). 
The Jeans mass decreases as the density increases up to the point at which the 
gas stops being approximately isothermal and begins to heat up adiabatically. 
From there on, the Jeans mass starts increasing again with density. The density 
at which this occurs defines the minimum Jeans mass that the simulation has to 
resolve and, thus, the minimum number of particles required. Present-day com- 
puter capabilities allow the use of several million SPH particles per simulation, a 
number that is nominally adequate for resolving the density at which adiabatic 
heating switches on. Using parallel computing facilities presently available, such 
computations take several months for each evolutionary model. 

The assumption of driven as opposed to decaying turbulence is important 
in these and similar simulations (cf. Heitsch, Mac Low & Klessen 2001, Li & 
Nakamura 2006). 

Indeed, Li & Nakamura (2006) using the best 3D MHD simulation tech- 
nique to date have shown that initial turbulence in a cluster forming region is 
quickly replaced by motions generated by protostellar outflows. This protostel- 
lar outflow-driven turbulence can keep a protocluster clump close to virial equi- 
librium long after the initial turbulence has decayed away. This may explain 
the early molecular cloud observations of Bertoldi & McKee (1992) and lend 
support to the recent equilibrium cluster formation models of Tan, Krumholz & 
McKee (2006) . It could even imply that the stellar IMF is regulated by outflow 
feedback (Silk 1995, Adams & Fatuzzo 1996). However, the ultimate test of 
these predictions has not yet been made and requires much more observational 
work. At the same time, competing, more violent cluster formation models have 
been proposed in the literature (e.g. Bonnell, Bate & Vine 2003), and the jury 
is still out. Essentially, the question is whether star formation is slow or fast 
(Ballesteros-Paredes, Hartmann & Vazquez-Semadeni 1999, Glover & Mac Low 
2007). 
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3.3 The Collapse Phase 



Gravity plays the dominant role in star formation. To form a star, gravity must 
overcome pressure, magnetic forces, internal turbulence, and rotation. In the 
simplest case of gravity versus gas pressure, one defines the Jeans mass. 



M 



Jeans ■ 



1.1 Mp 







-1 3/2 r 



Lgas 



10 K 



10-19 g cm-^ 



-1/2 



as the smallest mass for which gravity can become dominant. The normalization 
is consistent with typical initial conditions. 

Turbulence as a repulsive force will exceed gas pressure if motions are super- 
sonic. Unless continually replenished, however, supersonic turbulence dies out 
on a dynamical timescale (see, e.g.. Stone, Ostriker & Gammie 1998; Clark & 
Bonnell 2005; Kritsuk, Norman & Padoan 2006; for a summary of earlier work, 
see Mac Low & Klessen 2004). 

As Shu, Adams & Lizano (1987) point out, once gravity dominates pressure 
and magnetic forces in an optically thin gas capable of radiating compressional 
heat, it remains dominant and the gas collapses on a free-fall timescale 

tff~2.1xl05yr — / ^3 

[10 g cm 

That is, the densest parts collapse the fastest, and the Jeans mass decreases 
during collapse. The gas collapses nonhomologously until the densest parts be- 
come optically thick, allowing the gas to heat up adiabatically and to increase 
the gas pressure dramatically. Rotational (centrifugal) forces increase during 
gravitational collapse owing to conservation of angular momentum, so flattened 
structures and accretion disks (e.g., Black & Bodenheimer 1975; Terebey, Shu 
& Cassen 1984; Yorke & Bodenheimer 1999) are expected phenomena of gravi- 
tational collapse. 



3.4 The Accretion Phase 

3.4.1 Formation of the Hydrostatic Core 

The formation of low-mass stars is explained by the nonhomologous collapse 
of a slowly rotating fragment of molecular material; the collapse is stopped 
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in the central regions when the object becomes optically thick (Larson 1969, 
Woodward 1978, Winkler & Newman 1980, Masunaga, Miyama & Inutsuka 
1998; for a more recent review, sec Klein. Fisher & McKee 2004). There is 
a second, inside-out collapse in the center for Tgas ~ 2000 K when molecular 
hydrogen dissociates. When the second core is optically thick and thermally 
ionized, pressure forces are able to balance gravity on a dynamical timescale, 
and one speaks of an accreting quasi-hydrostatic core as more material rains 
down. When it forms, the second hydrostatic core contains somewhat more 
than a Jupiter mass and has a radius on the order of (3 - 5) (Winkler & 
Newman 1980, Masunaga & Inutsuka 2000). 

As long as material continues to flow onto the quasi-hydrostatic core, the 
core grows in mass. Simultaneously, it contracts on a thermal Kelvin-Helmholtz 
timescale tkh ^ GM^/R*L, toward hydrogen-burning densities and tempera- 
tures. Figure 4 shows the Kelvin-Helmholtz timescale to the ZAMS as a func- 
tion of stellar mass. Except for the growth of mass and angular momentum, 
the interior regions of the hydrostatic core are generally ignorant of and not 
influenced by the details of the accretion process. This allows us to separate 
the problem of star formation into several distinct parts: (a) the evolution of 
the central core, (&) the details of transporting material from the disk onto the 
core, (c) transporting material inward within the disk, and (rf) accretion onto 
the disk. Fundamental differences between low-mass and high-mass star forma- 
tion can be attributed to differences in the above processes and, in particular, 
to signiflcant differences in the timescales involved and in the local radiative 
environment. 

3.4.2 Evolution of Accreting Cores 

Accretion onto the (proto)star will affect the core's outer (atmosphere) regions 
and the overall spectral appearance. We estimate the thickness ARdyn of the 
outer stellar regions that are dynamically affected by accretion to be: 

ARdyn Ms-acc'S-dyn ts_dyn 



R, M* ts —acc 

where ts-dyn — R*/vesc, the stellar dynamical timescale, is the time for an 
acoustic wave to cross through the core or the orbital period of a body just above 
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the core's surface (generally less than a day). Because ts-dyn <<< ts-acc = 
M*/Ms-acc (the accretion timcscale), ARdyn/R* is extremely small (< 10~^ for 
accretion rates Mg-acc < 10^^ Mq yr^^). 

It is not easy to estimate the size of tlie region thermally affected by the 
accretion flow onto the (proto)star, because the radiation transfer in an optically 
thick plasma depends on details of the complex accretion geometry. Exactly how 
material flows from the disk onto the (proto)star remains an unsolved theoretical 
problem. 

One can expect a (proto)star's location in the Hertzsprung-Russell (HR) 
diagram to be strongly affected by accretion whenever the Kelvin-Helmholtz 
timescale tkh exceeds the accretion timescale Tacc- However, even when this 
criterion is not fulfilled, the Kelvin-Helmholtz timescale for the outer regions 
of a star can locally exceed the accretion timescale, and accretion can affect 
the radius, effective temperature, and luminosity of the star (Kippenhahn & 
Meyer-Hofmeister 1975). 

Kippenhahn & Meyer-Hofmeister (1975) consider the case of mass transfer 
within a close binary system. At high accretion rates Ms-acc ^ 10~^ Mq yr~^, 
accreting main-sequence stars of mass > 5 Mq bloat up to radii exceeding 10 
times their ZAMS values. Although not strictly applicable to the case of accret- 
ing premain-sequence stars, one can expect the same qualitative effect of high 
accretion rates in this mass range. These results demonstrate the importance 
of accretion for the appearance of the central hydrostatic cores. 

For low-mass stars, it is generally accepted that material is transported onto 
the central core through an accretion disk (e.g., Shu, Adams & Lizano 1987). 
The net gain per unit time of gravitational potential energy of the accreted 
material GM^Mg-acc/R-* is partly converted into rotational energy of the core 
and disk (-^1/4) and partially converted into heat (^3/4), which is radiated 
away (Yorke & Kriigel 1977). Part of the gravitational energy is converted into 
heat in a series of disk accretion shocks (Yorke & Bodenheimer 1999); part of it is 
converted into heat within the disk by the same viscous processes that transport 
angular momentum outward and allow the radial flow of material inward; and 
part of it is converted into heat in the accretion flow and shocks/relaxation 
zones on the (proto)stellar surface. 
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As mentioned above, the details of how the gas is ultimately transported 
from the disk onto the core are still unclear. For low-mass stars it has been 
postulated that magnetically focused flows (Blandford & Payne 1982) and/or 
accretion columns (Konigl 1991, Edwards ct al. 1993) and/or X-winds (Shu ct 
al. 1994, 1995) arc involved. Because 50% (~ 3/8 GM,Ms-acc/R*) of the total 
gravitational energy is converted into heat and radiated away within 1 of 
the core, one can expect higher temperatures and luminosities as either or 
Ms-acc increases or as the core contracts. The temperature of this hot shocked 
gas 



10^ K 
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10 Rr; 
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is sufficiently high to produce X-rays even for low-mass stars. 

The above argument is only slightly modified when the core and disk produce 
a wind or outflow. We then speak of a net mass accretion of Ms_acc — Mg-wind- 
A portion of the core's rotational energy and angular momentum can be con- 
verted into kinetic energy and angular momentum of the wind. 

A necessary condition to accrete sufficient material to produce a massive 
star is thus, 

■ dt'>8MQ. (1) 
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Ms_acc(t') - Ms_^i„d(t') 



That is, the accretion rate Ms-acc onto an embryo object must exceed the 
outflow rate Ms-wind during a significant proportion of the formation process. 
For this to occur, the acceleration owing to gravity must exceed the outward 
directed radiative acceleration of the accreting core. Whereas gravity GM*/r^ at 
each radial point in an envelope increases linearly with core mass, the radiative 
acceleration of dusty material KL/47rr^c is proportional to the core's luminosity. 
This increases as a high power of stellar mass. A lower limit to the core's 
luminosity is the ZAMS luminosity (see Figure 4). 

Thus, to allow infall, we require as a necessary condition 



KesL/Am'^c < GM*/r^ 
with L = L» + Lace 5 which translates into 

Keff < 130cm^g~^ 
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where the effective opacity for radiative acceleration of accretable material is 
defined as 



here k^J^ is the frequency-dependent gram-opacity (cgs-units: cm^g~^) of the 
material subject to radiative acceleration and F = Jg°° F^, du is the radiative 
flux. (For a more complete discussion of the definition of k]^, including effects 
of nonisotropic scattering and photoejection of particles from the dust, see Yorke 
1988.) The (proto)star's luminosity is given by the sum of its intrinsic luminosity 
L» and the luminosity Lace emitted by the dissipation of kinetic energy of the 
material being accreted. 

For dusty gas, K^f^ is strongly frequency-dependent. Depending on the hard- 
ness of radiation emitted by the central source and the accretion shocks, dusty 
gas could be repelled from the star. The ISM extinction in the J-band, for 
instance, corresponds to k = 130cm^ g~^. In Figure 5 we indicate values of 
Kes for dusty gas at solar abundances {n^f''^ taken from the Preibisch et al. 
1993 dust model with and without ice-coated grains) illuminated by 600 K and 
6000 K black bodies [i.e., Fi, cx B^(T*)], respectively, and for a fully ionized 
dustless hydrogen plasma. 

Figure 6 displays the mean effective opacity of dusty gas as defined by 
Equation 4, using the Preibisch et al. (1993) dust model with and without 
ice-coated grains {blue and red curves, respectively) and assuming blackbody 
radiation at the temperature Trad- Using the mass scale on the right of Figure 6 
we note that the net force on dusty gas surrounding a deeply embedded 100 Mq 
main-sequence star is directed toward the star as long as the star appears to 
the dust to be a <, 50-K source. By contrast, for an unobscured 3M0 main- 
sequence star with T^s ^ 10'' K the net force on dusty gas would be directed 
away from the star. 

3.5 Overcoming Radiative Acceleration 

Next we discuss how Equation 2 can be satisfied, i.e., accretional growth of 
an already existing stellar embryo can be enabled. At least one of the following 
conditions must be met: (a) k^s must be sufficiently low, i.e., significantly lower 
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than its ISM value for optical/UV radiation; (6) the total luminosity L must be 
reduced; or (c) gravity (i.e., the stellar mass M* ) must be increased. Below, we 
discuss each of these three possibilities. 

3.5.1 Reduce Keff 

Owing to the strong frequency dependence of dust opacity, Keff can be signif- 
icantly lower than its ISM value if the radiation field seen by the accreting 
material is that of a cold, embedded object. When embedded the protostel- 
lar radiation field is shifted from the optical/UV - where the dust absorbs the 
photons - into the far-IR where the dust reemits the absorbed energy. Alter- 
natively, Keff can also be reduced if the average size of dust grains increases 
(but remains compact rather than becoming fractal) or if most of the dust 
is destroyed. In their pioneering efforts, Kahn (1974) and Wolfire & Cassinelli 
(1987) studied the ID, spherically symmetric accretion problem for massive star 
formation with an emphasis on the dust opacity. Kahn concluded that a 40- 
Mq star could be formed by spherically symmetric accretion, but his assumed 
dust destruction temperature Tgub = 3600 K is too high by a factor of about 
two (implying for the outer regions greater protection from radiation pressure) . 
Furthermore, his assumed grain opacities, which can be written in the form 
K ~ 100cm^g~^(Tdust/9000K)^, were somewhat too low. 

In a more careful treatment of the dust (but still assuming steady-state 
spherically symmetric infall), Wolfire & Cassinelli (1987) concluded that very 
massive stars can form only if the dust has been significantly modified. By 
contrast, Yorke & Kriigel (1977) showed in a hydrodynamical simulation that 
spherically symmetric accretion must be nonsteady for the high-mass case. They 
were able to produce stars of masses 17 Mq and 36 Mq from clouds of masses 
50 Mq and 150 Mq, respectively, in a highly variable accretion flow. 

These early attempts to explain high-mass star formation suffer significantly 
from the fact that spherically symmetric infall was assumed. 

The effective opacity of the accreting material can also be reduced by den- 
sity inhomogeneities resulting from the photon bubble instability. The radiation 
escapes readily through the gaps between the shocks that are driven by distur- 
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bances in the radiation flux (Turner, Quataert & Yorke 2007). Another possibile 

explanation for reducing the effective opacity is the accretion of optically thick 
blobs or fingers. These can be expected to form via Ray Icigh- Taylor instabilities 
in the radiation-inflated cavities produced by luminous (proto) stars (Krumholz, 
McKee & Klein 2005a; Yorke 2002). In this case, 

KeS = TrRyob/^blob- 

As a particular subset of this family of solutions, Bonnell, Bate & Zinnecker 
(1998) considered building up massive stars by coalescence of intermediate-mass 
stars within a very dense protostellar cluster (see also Bonnell & Bate 2002, 
Zinnecker & Bate 2002). We return to this issue in more detail below. 

Opacity modifications owing to coagulation of dust and dust destruction 
processes during the collapse phase were calculated by Suttner & Yorke (2001) 
for three different detailed dust models: compact spherical particles, fractal 
BPCA (ballistic particle-cluster agglomeration) grains, and fractal BCCA (bal- 
lastic cluster-cluster agglomeration) grains. (BPCA dust consists of grains with 
widely different grains sizes, whereas BCCA dust is formed by coagulation of 
similar-sized grains). Assuming axial symmetry, Suttner & Yorke followed the 
dynamics of gas and 30 individual dust components. They found that even 
during the early collapse and the first ^ 10^ yr of dynamical disk evolution, the 
initial dust size distribution is strongly modified. Close to the disk's midplane, 
coagulation produces dust particles of several tens of /im in size (for compact 
spherical grains) up to several millimeters in size (for fluffy BCCA grains). In 
contrast, in the vicinity of the accretion shock front located several density-scale 
heights above the disk, large velocity differences inhibit coagulation. Dust par- 
ticles larger than about 1 /im segregate from smaller grains behind the accretion 
shock. Owing to the combined effects of coagulation and grain segregation the 
IR dust emission is modified. Within the accretion disk, an interstellar medium 
dust size distribution (Mathis, Rumpl & Nordsieck 1977) provides a poor de- 
scription of the general dust properties. Nevertheless, the radiative force acting 
on material infalling from the envelope to the disk is hardly affected by coagu- 
lation (see also Weidenschilling & Ruzmaikina 1994). 
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3.5.2 Reduce the Effective Luminosity 

Nakano, Hasegawa & Norman (1995) and Jijina & Adams (1996) pointed out 
that, because we expect accretion to proceed through an accretion disk, radia- 
tion pressure could blow away the tenuous polar regions but not the massive disk. 
Yorke & Bodenheimer (1999) and Yorke & Sonnhalter (2002) studied this ef- 
fect quantitatively and substantiated this claim through numerical simulations. 
They found that, whereas the central object may emit radiation isotropically, 
the radiation field quickly becomes anisotropic farther from the center. For an 
outside observer, and in particular, for a dust grain attempting to accrete onto 
an existing protostellar disk, the radiative flux close to the equatorial plane can 
be much smaller than the component parallel to the rotation axis. This so-called 
flashlight effect (beaming of radiation in the polar direction) occurs whenever a 
circumstellar disk forms. 

Yorke & Bodenheimer (1999) considered the evolution of a lO-M© molecular 
fragment that produced a 8.2-M0 star (their case F). Figure 7 depicts an 
intermediate stage of the evolution, after 7 of material have accreted onto 
the central protostar, while 2.8 M0 have accumulated in a circumstellar disk and 
0.2 Mq still reside in an larger extended inf ailing envelope. Inner and outer disk 
accretion shocks (density and temperature discontinuities) are visible above and 
below the disk. The dark red regions with an opening angle of about 45° show 
preferential heating of polar regions by the protostar. The Ught red regions 
within 45° of the equator beyond ~ 1000 AU indicate shadowing by the disk. 

Yorke & Bodenheimer (1999) estimated that the edge-on and pole-on bolo- 
metric fluxes can differ by more than a factor of 30 after about one-half of the 
mass of a 2-M0 collapsing protostellar clump has accreted onto the protostar. 
The difference in radiative acceleration is much greater than this factor of ^ 30, 
however, because the edge-on flux is dominated by the far-IR, which is far less 
effective at radiatively accelerating dusty gas than mid- and near-IR light seen 
pole-on. Yorke & Sonnhalter (2002) showed by calculating frequency-dependent 
radiation transfer that the flashlight effect is further enhanced by the fact that 
the central star's optical and UV radiation blows out material in the polar di- 
rection, reducing back-scattering of radiation toward the disk. A similar effect 
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has been reported by Krumholz, McKee & Klein (2005a) from polar cavities 
blown free by outflows originating close to the star. Thus, photons emitted or 
scattered into these directions will not hinder accretion of material within the 
disk or material in the disk's shadow regions. 

Although the flashlight effect allows dusty material to come close to the 
central source via a circumstellar disk, the material to be accreted eventually 
encounters optical and UV radiation from the central source. For this material 
to be accreted rather than blown out by radiation, the dust must be largely 
destroyed or it must have coagulated into larger particles so that the opacity is 
dominated by the gaseous component. 

Even though no massive disk has yet been directly observed around a main- 
sequence O star, there is much indirect evidence that such disks exist (see the 
review by Zhang 2005 and our more detailed discussion in Section 6). A 
compelling argument that disks exist during the early phases of massive star 
formation is the observation of massive bipolar outflows. Such massive outflows 
are probably powered by disk accretion, and, as in their low-mass counterparts, 
the flow energetics appear to scale with the luminosity of the source (see Shep- 
herd & Churchwell 1996a,b; Richer et al. 2000; Henning et al. 2000; Wu et al. 
2005). 

If the primary source of the massive star's material is accretion from the 
surrounding molecular core, then a circumstellax disk should be the natural 
consequence of the star-formation process even in the high-mass case. However, 
it should be diflicult to observe disks around massive stars. The high far UV and 
extreme UV fluxes associated with high-mass stars will begin to photoevaporate 
the disks on timescales of ~ 10^ yr (Hollenbach, Yorke & Johnstone 2000) . The 
results will be observable as deeply embedded ultracompact Hll-regions with 
comparable lifetimes (Richling & Yorke 1997). The fact that disks around O 
stars photoevaporate so quickly provides negative feedback for disk accretion. 
This limits the build-up of more massive stellar objects and may even imply an 
upper mass limit for star formation. 
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3.5.3 Increase Gravity 

For completeness we mention the fact that the gravitational acceleration is en- 
hanced with respect to radiative acceleration when massive stars form within a 
dense cluster of not-so-brightly-radiating objects. In this scenario, one requires a 
density-peaked cluster of low-mass objects embedded within a molecular cloud, 
with Pobjects ^ Pgas- The effective gravity near the cluster's center is enhanced 
relative to an isolated molecular cloud without the cluster and relative to off- 
center regions of the molecular cloud (cf. Keto 2002) . If this were the only way 
to form massive stars, isolated massive stars would exist only in very exceptional 
cases. 

3.6 Stelleir and Protostelleir Luminosity Evolution 

Because the luminosity is so critical during accretion up to high stellar masses, 
one must also consider the luminosity evolution of the accreting object. As 
discussed by Behrend & Maeder (2001) and Yorke (2002), Figure 8 shows 
that protostars do not evolve along premain-sequence tracks until they land on 
and remain at a unique spot on the main sequence where hydrogen burning 
starts. Rather, they reach the main sequence - that is, the hydrogen-burning 
state - well before they have finished accreting mass. After that, they con- 
tinue growing in mass and evolve up the main sequence until they run out 
of material to accrete. This means that an initially low-mass object that gains 
mass through accretion evolves substantially differently in the HR diagram than 
would a nonaccreting premain-sequence star of the same final mass. 

Figure 8 shows a number of protostellar evolutionary tracks in the HR di- 
agram, assuming a sequence of conceivable accretion rates increasing by factors 
of 10, but each constant in time. These accreting tracks were calculated by 
Yorke (2002) and are qualitatively similar to the more detailed stellar evolution 
calculations by Norberg & Maeder (2000), Behrend & Maeder (2001), and - for 
the lower accretion rates - Palla & Stabler (1992). Differences to the two former 
investigations can be attributed to Yorke's (2002) starting mass, 0.1 M© instead 
of 1 Mq, and the differing accretion rates. In all cases published to date, not only 
do the tracks of accreting objects consistently lie slightly below the equilibrium 
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deuterium-burning birthline, but the qualitative effect of more rapid accretion is 

to shift the tracks away from the birthUne. (In this context, the birthline is de- 
fined as the point on the Hayashi track at which equilibrium deuterium burning 
begins in a nonaccreting star. Palla & Stahlcr (1993), however, use the concept 
of birthline as the evolutionary track of accreting prcmain-scqucncc stars in the 
HR diagram. Deuterium burning of newly accreted material keeps these tracks 
well above the main sequence until accretion stops.) The tracks of accreting 
stars eventually converge to the main sequence and follow along the ZAMS as 
more material is added. For example, at an accretion rate of 10~'^MQyr^^, 
hydrogen burning begins at t ~ 1.3 x 10^ yr, after ~ 13 Mq have been accreted. 

The tracks discussed above do not display the degree of bloating seen by Kip- 
penhahn & Meyer-Hofmeister (1975) for accreting main-sequence stars M* ~ 5- 
10 Mq. Yorke's (2002) simplifying assumption of thermally adjusted (pre)stellar 
objects is not strictly fulfilled, especially for the Mg-acc = 10~^ yr~^ case. 
Norberg & Maeder (2000) did not consider accretion rates Mg-acc > 10~* Mq yr~^ 
and Behrend & Maeder (2001) attained high accretion rates Ms-acc ^ 10~^ Mq yr~^ 
for high-mass objects > 25 only. Palla & Stabler (1992, 1993) found sig- 
nificant bloating for the highest accretion rate lO~''M0yr^^ they considered, 
attributable to shell deuterium burning. We expect that (proto)stars accret- 
ing at extremely high rates are not fully thermally adjusted and thus bloat up. 
The degree of bloating depends on the accretion rate and mass of the accreting 
object (Kippenhahn & Meyer-Hofmeister 1975). 

We remind the reader that these theoretical tracks in the HR diagram do 
not refiect the actual observable bolometric luminosities of accreting protostaxs 
nor do they reflect the observable effective temperatures. Plotted in the HR 
diagram are the intrinsic (proto)stellar properties, from which the observed 
properties must be derived by multidimensional radiative transfer calculations 
(e.g. Indebetouw et al. 2006). Note that much of the accretion luminosity, 
Lacc5 is indistinguishable from the intrinsic luminosity, L*, of the star. Also, 
owing to the existence of shock fronts and their postshock relaxation zones, 
strong deviations from a stellar spectral energy distribution can be expected. 
The importance of accretion luminosity. Lace ~ GM»M/R», is seen in Figure 5. 
At high accretion rates the L/M ratio can attain high values even for low- mass 
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cores, so that radiative acceleration can clear out part of the infalling material 
in the polar regions. 

What order of magnitude of mass accretion rate can be expected? In order 
to produce a star of mass within, say 200,000 yr, an average accretion rate 
of 5 X 10~^ M0 yr~-'^[M*/M0] is necessary. Assuming this average accretion 
rate, wc note that during the main accretion phase the himinosity of low-mass 
stars is dominated by accretion luminosity, whereas for high-mass stars the 
luminosity is initially determined by accretion but is eventually dominated by 
the intrinsic stellar luminosity. Of course, the actual accretion rate may vary 
strongly from this average value. The maximum sub-Eddington accretion rate 
possible onto a core hydrogen-burning star, assuming electron scattering and the 
effects of both the intrinsic stellar luminosity and accretion luminosity, can be 
inferred from Figure 5. Whereas the dotted curve corresponding to an accretion 
rate of lO~^M0yr~^ lies below the value permitted by electron scattering, an 
accretion rate ten times higher would clearly lie above this value everywhere. 
Thus, accretion of ionized material onto a stellar core at a rate ^ 10~^ Mq yr~^ 
implies super-Eddington accretion. 

3.7 Stelleir Evolution Beyond the Zero Age Main Sequence 

As hydrogen burning proceeds, the stars begin to evolve away from the ZAMS 
even as they accrete material through a disk and lose material through a stellar 
wind. The star's evolution strongly depends on its mass loss rate and on the 
internal mixing that is induced by rotation. These effects must be taken into 
account when modeling these stars and their effect on their environment. 

If the stars attain their final mass and become optically visible on a timescale 
much shorter than their main-sequence lifetime, they will not have evolved far 
from the ZAMS. If, however, the stars first become optically visible after signifi- 
cant stellar evolution, the traditional concept of the ZAMS as the starting point 
for newly formed massive stars is fiawed. To illustrate this point we compare 
in Figure 9 the current best estimate of the locations in the HR diagram of 
spectral types 03 through 09.5 for luminosity classes V, III, and I with theoret- 
ical isochrones and evolutionary tracks. It is not surprising that the O dwarfs 
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do not fall on the ZAMS, because we are dealing with average properties over 
a range of ages. A quantitative explanation of the magnitude and systematic 
variation of the offset for early to late O dwarfs, for example a shift of average 
ages from ^ 1 to 5Myr, is still lacking (Martins, Schacrer & Hillicr 2005). Such 
a shift could occur, because stars of later spectral types have weaker winds and 
lower UV fluxes and therefore remain embedded in their parental cloud longer. 

Table 3: 0-stars spectral type versus mass for different luminosity classes (Mar- 
tins, Schaerer & Hillier 2005) 



SpT 




Luminosity clciss 






V 


III 


I 


BO 








09.5 


16.5 


21.0 


30.4 


09 


18.0 


23.1 


32.0 


08 


22.0 


26.9 


36.8 


07 


26.5 


31.2 


40.9 


06.5 


29.0 


33.7 


43.1 


06 


31.7 


36.4 


45.8 


05 


37.3 


41.5 


50.9 


04 


46.2 


48.8 


58.0 


03 


58.3 


58.6 


66.9 


02.5 









Note: masses in solar units 



We find it useful to end this stellar evolution section with a compilation 
(Table 3) of the masses of massive stars as a function of spectral type (from 
09.5 through 03) and luminosity class (V, III, and I). 

3.8 Ionization Evolution and Cloud Disruption 

Because massive stars are built by accretion, whether from a monolithic collapse 
or competitively, they will produce an increasing number of hydrogen ionizing 
photons even while they grow as they become hotter and more luminous. Thus, 
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in spite of accretion, hydrogen-ionizing and helium-ionizing fluxes similar to 
those expected from ZAMS stars are likely (sec Figure 10). Indeed, the pro- 
cess of accretion itself is likely to produce hard ionizing radiation in a scries of 
accretion shocks close to the stellar surface. Moreover, powerful winds interact- 
ing with surrounding material will produce strong shocks and hard radiation. 

Thus, material in the immediate vicinity of the accreting massive star, in 
particular material in the circumstellar disk and in any nearby star-disk systems, 
can be ionized. In the case of circumstellar disks, a thin layer of ionized gas 
on the surface of the disk results. Close to the star, where the escape velocity 
is much greater than the sound speed, ionized gas remains bound. Ionized gas 
outside a radius revap, that is, the radius where the sound speed exceeds the 
escape velocity, expands outward. The situation depicted in Figure 11 is likely 
to result. 

The photoionization of circumstellar disks produces an outflowing disk wind 
of ionized gas that interacts with the stellar wind. Depending on its strength 
(and the photoionization rate) , the stellar wind will be collimated by the ionized 
disk wind (Richling & Yorke 1997). The stronger the stellar wind (or the weaker 
the ionizing photon flux and hence the ionised disk wind), the less collimated 
the stellar wind will be. 

This is different from the case of collimated jets and outflows due to X-winds 
and disk winds in young low-mass stars, where magnetic fields and forces are 
involved (Shu et al. 1995). In the high- mass case, although magneto-centrifugal 
forces cannot be dismissed from influencing the collimation process, we can 
have rather collimated outflows without magnetic flelds. This being said, we 
suspect that magnetic fields do play a role not only in low-mass but also in 
high-mass star formation, not least because there are by now at least two cases 
where strong kilo-Gauss stellar magnetic fields have been inferred from Zeeman 
spectro-polarimetry, including the oblique magnetic rotator Ori C (Stahl 
et al. 1996, Donati et al. 2002) and the Of?p star HD 191612 (Walborn et al. 
2003, Donati et al. 2006). The mismatch in 0^ Ori C between the magnetic field 
direction and the stellar spin axis could be an indication that this well-known 
massive star formed in a coUisional process from a merger of two lower mass 
stars. Finally, we mention the detection of hard X-ray emission from the W3 
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massive star-forming region, an indication that embedded massive stars must 
be magnetically active (Ho&ier et al. 2002). 

Hoare et al. (2007) examined the suggestion first made by HoUenbach et 
al. (1994) and later calculated numerically by Yorke & Welz (1996) and Rich- 
ling & Yorke (1997) that photoevaporating disks can explain the existence of 
unresolved ultraconipact HII regions (sizes ^O.lpc). Many of these have since 
been resolved with long-baseline radio interferometry and have been reclassifed 
as hypercompact HII regions (sizes ~0.01pc; see Kurtz 2005). The latter most 
likely represent the individual photoevaporating disks associated with individ- 
ual deeply embedded O stars (Keto 2007). The phase of photoevaporating disks 
is the first evolutionary stage after the ionizing flux has turned on. It remains to 
be seen whether the thermally evaporating flow interacting with a stellar wind 
can produce the high velocities seen in the broad H92a recombination line seen 
toward hypercompact HII regions (Sewilo et al. 2004). 

Hoare et al. (2007) concluded that in a later evolutionary phase most ultra- 
compact HII regions need to be interpreted as external photoevaporation of a 
molecular clump. The longevity of these HII regions is explained by a relative 
motion of the ionizing source and its stellar wind into the clump, producing a 
combination of bow shock and unipolar flow (a so-called champagne flow, see 
the review by Yorke 1986). If this is the case, then the transition from hyper- 
compact to ultracompact HII regions may mark the end of the existence of the 
accretion disk. The ultracompact HII region phase then marks the beginning of 
the dissipation of the molecular cloud as the stellar wind and ionizing radiation 
of the newly formed massive star are able to interact with lower density, more 
poorly gravitationally bound molecular material. 

An important point to make in the evolution of an ultracompact HII re- 
gion is that the gravitational force of the star(s) responsible for the HII region 
should be included (Keto 2002). The gravitational attraction of the star(s) can 
maintain the accretion flow within the ionized gas and prevent the HII region 
from expanding hydrodynamically. This is true as long as the radius of the 
ionization equilibrium is smaller than the radius where the sound speed of the 
ionized gas, about lOkms"^, approximates the escape velocity (concept of the 
trapped HII region). Indeed, observations of the H66ci! recombination line from 
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the ultracompact HII region GlO.6-0.4 powered by a compact and very lumi- 
nous (~ 10® L0) cluster of newly formed massive stars surprisingly show inward 
motion (Keto 2002). This implies that, despite the high luminosity and ioniz- 
ing radiation of several O stars, neither radiation pressure nor thermal pressure 
has reversed the accretion flow (Keto & Wood 2006). The ram pressure of the 
observed accretion flow, with a rate of M* = 10""^ M0 yr~^, can conflne the bub- 
bling HII region and allow the massive stars to continue growing in mass - from 
ionized gas! One may speculate whether the ionization surrounding the moder- 
ately massive stars helps rather than hinders their further growth by accretion 
(Keto & Wood 2006), the reason being that ionized gas can couple strongly with 
any magnetic fields. Magnetic fields can transport angular momentum outward, 
thereby allowing a high rate of infall to smaller radii to be maintained. Eventu- 
ally, of course, with the stellar masses steadily increasing, the ionizing flux will 
increase enough (roughly ocM^, see Figure 10) so that the radius of ionization 
equilibrium grows beyond the critical radius within which the gas is gravita- 
tional bound. At this point the HII region will burst free and start to disrupt 
the dense cloud, preferentially unidirectionally in a champagne flow, revealing 
the initially deeply embedded cluster (both high- mass and low-mass stars). 
Interim Summary 

Despite the obstructing effects of radiation pressure on dust, which would pre- 
vent massive stars with masses of more than 30 - 40 Mq to be formed in spher- 
ically symmetric collapse and accretion models, these and more massive stars 
can be formed via accretion through a circumstellar disk (i.e., in 2D and 3D 
models). An accreting star quickly evolves to the main sequence after about 
10 M0 have been accreted, but the star is not yet observable. Still obscured by 
the material in its vicinity, its appearance will be that of a hypercompact HII 
region. Radiative acceleration, photoevaporation, and stellar winds eventually 
destroy the accretion disk, but prior to this, accretion onto the star provides 
an additional source of luminosity. (Realistically, this accretion is expected to 
be highly variable and episodic.) A powerful radiation- or magnetically driven 
outflow in the polar directions and a puffed-up (thick) disk result from the high 
luminosity of the central source. The details of how disk material ultimately 
flows onto the star are still unclear - as is often the case for accretion disks. 
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4 MASSIVE STAR FORMATION: COMPET- 
ING CONCEPTS AND CALCULATIONS 

Three different concepts describing the origin of massive stars have been dis- 
cussed in the recent literature, each of which may occur in nature, depending 
on the initial and environmental conditions for the parent molecular clouds 
(e.g., the Mach number of supersonic turbulence in the clouds and the external 
pressure in the ISM). These arc (1) monolithic collapse and disk accretion, (2) 
competitive accretion and runaway growth, and (3) stellar collisions and merg- 
ers. Whether the parent molecular cloud clumps are magnetically subcritical or 
supercritical (Shu, Adams & Lizano 1987, Crutcher & Troland 2007) may play a 
crucial role in determining whether massive stars form (1) in isolation or (2) are 
strongly clustered. Furthermore, in extreme star-forming environments (such as 
in massive protoglobular cluster clouds), the initial gas densities may have been 
so high that (3) stellar collisions become an unescapable ingredient in massive 
star formation. We discuss all these routes toward massive star formation in this 
section. In addition, we briefly discuss competing, but not mutually exclusive 
scenarios related to massive star formation stimulated by the pressure of ex- 
panding HII regions versus a supernova blast wave. Similarly, we highlight the 
two competing, hut not mutually exclusive processes accounting for runaway 
OB stars. Finally, the important concept of mass segregation (the fact that 
massive stars are often concentrated in the centers of star clusters) is reviewed 
in terms of nature (a birthmark) or nurture (subsequent dynamical evolution). 

4.1 Monolithic Collapse and Disk Accretion 

Yorke & Sonnhalter (2002) consider the collapse of isolated, rotating, nonmag- 
netic, massive molecular cores of masses 30 Mq, 60 Mq, and 120 Mq using a 
frequency-dependent radiation hydrodynamics code. The flashlight effect dis- 
cussed in Section 3.5 allows material to enter into the central regions through 
a disk. For massive stars, it is important to take into account the frequency- 
dependent nature of the opacity and the flux within the disk rather than as- 
suming either Rosseland or Planck gray opacities. For their 60 Mq case, Yorke 
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& Sonnhalter find that 33.6 are accreted in the central regions as opposed 
to 20.7 Mq in a comparison gray calculation. Because these simulations cannot 
spatially resolve the innermost regions of the molecular core, they cannot distin- 
guish between the formation of a dense central cluster, a multiple-star system, 
or a single massive object. They also cannot exclude significant mass loss from 
the central object(s) that may interact with the inflow into the central grid cell. 
With the basic assumption that all material in the innermost grid cell accretes 
onto a single object, they are only able to provide an upper limit to the mass 
of stars that could possibly be formed for the cases considered. 

Note that the ~ 43 star formed during the collapse of the 120 molec- 
ular clump (Yorke & Sonnhalter 2002) does not represent an upper limit to 
the stellar mass that is enabled by the flashlight effect. Larger initial masses 
or a more focused flow along a fllament could conceivably lead to even more 
massive stars. In an ongoing simulation of the collapse of 100- and 2OO-M0 
clouds, Krumholz, Klein & McKee (2005) report that ^ 27M0 25 for 
the turbulent case) accreted onto a stellar core. As in the simulations of Yorke 
& Sonnhalter (2002), a disk formed around the accreting (proto)stars. 

One can speculate on the effect outflows have on the accretion through an 
accretion disk. The inner part of the accretion disk could well look like the 
configuration shown in Figure 11. Radiation and the stellar wind from the 
central star (presumably already hydrogen burning) evacuate a cavity in the 
polar direction. At the interface between the supersonic outflowing stellar wind 
and the denser subsonic HII disk atmosphere, some disk material will be re- 
moved, but this cannot prevent inward flow of disk material. Inward radial flow 
of dusty molecular gas is allowed in the equatorial plane of the disk as angular 
momentum is transfered outward. 

In cases of magnetized disks with high radiation densities, photon bubbles 
can lead to more eflicient vertical transport of radiation in the disk (Gammie 
1998; Begelman 2001; Blaes & Socrates 2003; Turner et al. 2005; Turner, 
Quataert & Yorke 2007). This can explain super-Eddington accretion in a vari- 
ety of luminous systems, including accreting compact objects and very massive 
stars. 

Eflicient angular momentum transfer could result from weak magnetic flelds 
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in the disk (Balbus & Hawley 1991, Hawley & Balbus 1991, Balbus 2003), 
from turbulence and/or spiral density waves (Bodenheimer 1995) excited by 
gravitational instabilities in nonmagnetized (Laughlin & Bodenheimer 1994) or 
magnetized (Fromang et al. 2004) disks, or from the tidal effects of nearby 
stars (Terquem 2001). Indeed, rapid accretion through a disk may be a direct 
consequence of having nearby companions. This may explain why massive stars 
are generally members of multiple systems. 

Once the disk material crosses rdust ~ 25 AU [M*/30 Mo]^-^, the radius of 
dust destruction, its opacity decreases and it is not easily stopped by radiation. 
It is, however, still unclear how the disk material ultimately flows onto the star. 
Surely, the disk puffs up close to the star, in analogy to the accretion disks 
in active galactic nuclei. Beyond rgvap ^ 13OAU[M*/3OM0], where the escape 
velocity is less than lOkms"^, the disk loses material via photoevaporation on 
a timescale of 10^ yr. This is of the same order as the accretion timescale. 
These competing effects (accretion and photoevaporation) will determine the 
final mass of the star and perhaps even the upper mass limit. 

4.2 Competitive Accretion and Runaway Growth 

Bonnell et al. (1997, 2001a) present the first 3D numerical simulations of the 
growth of stellar masses by competitive accretion in small young star clusters. 
The best way to visualize the idea of competitive accretion is to compare it 
to an economic model based on two complementary concepts; the real estate 
concept - "location, location, location" - and the capitalistic concept - "the 
rich get richer." The former reminds us that environmental infiuences can be 
very important, such that being in a fortunate location can significantly pro- 
mote growth, whereas the latter simply means that the gravitational attraction 
increases with success - in this case, with the increasing mass of the star. A 
protostar's ability to grow depends on the size of its accretion domain, i.e., the 
region from which gas can be gathered. A location in the center of a protostellar 
cluster is beneficial, as gas fiowing down to the center of the cluster increases 
the gas reservoir available to an individual star. The early birth of a protostar 
may also give it an unfair advantage in the competition to end up large and 
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massive. 

Imagine a large, dense molecular gas cloud with a number of protostel- 
lar seeds distributed inside the cloud that have initially condensed from some 
denser portions of the cloud. These condensations or cores subsequently have the 
chance to grow in mass by accumulation (accretion) of lower density cloud gas 
from their individual accretion domains (Larson 1978, Zinnecker 1982). These 
accretion domains are systematically larger for higher mass seeds. With increas- 
ing mass, the gravitational spheres of influence keep growing. 

Equally important, the amount of material that enters an individual accre- 
tion domain depends on the external environment. In isolation, it would simply 
be proportional to the mean gas density of the region. By contrast, in a larger 
scale potential such as that provided by a protocluster cloud, the gas density 
can become significantly larger in the center of the protocluster, as the gas set- 
tles into the deepest part of the potential, there to be accreted by the growing 
protomassive star. 

The accretion domain of an off-center protostar is tidally limited by the total 
mass in the inner part of the cloud, whereas the accretion domain of a protostar 

in the center of the cloud is the whole cloTid. Because the gas reservoir is limited 
(the cloud has a finite mass), the protostellar masses will eventually compete for 
cloud gas, especially after the accretion domains start to overlap. The action 
of the cluster to gather matter from larger distances and focus it toward the 
accreting stars, combined with the increasing accretion radii of these stars owing 
to their increasing mass, is what makes competitive accretion such a powerful 
mechanism. 

This is also true in a molecular cloud with hierarchical substructure, in 
particular for subcluster clumps of gas in a bigger protocluster cloud. The 
implication is that each subcluster clump is likely to have one most massive 
protostellar object in its center surrounded by a hierarchy of lower mass objects 
(Bonnell, Bate & Vine 2003). When and if those subcluster clumps merge to 
one big protocluster (see Figure 12), the complex accretion history of each 
protostar is no longer related to a single local cloud core (Schmeja & Klessen 
2004). In this scenario, the massive stars had accretion histories that were 
priviledged at every stage - from Jeans-instability protostellar birth to location 
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and density in the gas cloud, all factors affecting accretion were more favorable 
than the average. This is why massive stars are rare. The rarest and most 
massive of them probably formed in the most favorable conditions by runaway 
accretion until their gas reservoir was exhausted or dissipated (e.g., by ionization 
feedback, sec von Hocrncr 1968 and recently Clarke, Edgar & Dale 2005). 

The above model of competitive accretion has been critizcd by Krumholz, 
McKce & Klein (2005b) on the grounds that Bonnell's SPH simulations start 
from very strongly gravitationally bound protocluster clouds, while observa- 
tionally such clouds appear to be supported by turbulent motions. In other 
words, the simulations use a virial parameter Q! = Eturb/Egrav ^ 1 while molec- 
ular observations suggest a~ 1. Krumholz, McKee & Klein (2005b) argue from 
analytical considerations that protostellar masses cannot grow in such a turbu- 
lent medium, not even by a factor of two. The Bondi-Hoyle accretion rate is far 
too low owing to the high relative velocities between the accreting stars and the 
turbulent gas. In addition, radiative feedback from the incipient massive star 
may prevent Bondi-Hoyle accretion altogether (Edgar & Clarke 2004). 

In response to this criticism, Bonnell & Bate (2006) argue that global initial 
collapse versus quasi-cquilibrium support is not the issue, as their more recent 
SPH calculations start with initial conditions where the turbulent kinetic energy 
is close to the gravitational energy. In fact, the scaling laws of supersonic 
turbulence (Larson 1981) imply small turbulent velocity differences between 
protostars and their neighboring gas. This then allows for significant growth in 
stellar mass while more distant high velocity turbulent gas cannot be accreted. 
In other words, protostars are swept along with their neigboring gas for some 
time in a similar global motion before they hit more distant and thus less- 
correlated gas. The above Bondi-Hoyle accretion problem is also not as serious 
as suspected, considering that the most massive star forms in the center of the 
clump with little motion relative to the surrounding gas. 

Future detailed observations of the motions in massive (-^IOOOMq) star- 
forming clumps will discriminate between these different dynamical scenarios 
(small versus large gas motions relative to the protostars, turbulent cloud sup- 
port, or overall cloud collapse). A first step in this direction was made by 
Peretto, Andre & Belloche (2006) in their dust continuum and molecular study 
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of the NGC 2264 clumps, with the conclusion that the observations are most 
consistent with "a picture of massive star formation intermediate between the 
scenario of stellar mergers of Bonnell et al. (1998) and the massive turbulent 
core model of McKee & Tan (2003), whereby a turbulent, massive ultra-dense 
core is formed by the gravitational merger of two or more Class protostellar 
cores at the center of a collapsing protocluster" . 

4.3 Stellar Collisions and Mergers 

Historically, the original reason for proposing stellar collisions as a formation 
process for massive stars was twofold, (a) At the time, radiation pressure on 
dust was considered a severe hindrance to gas accretion (then assumed to oc- 
cur spherically symmetrically). Today this concern has gone away (see Sec- 
tion 4.1). (6) The packing of massive stars in dense clusters was too tight, 
and so there was concern that a sufficiently large gas reservoir for monolithic 
collapse was not available. This is still a concern, although not all massive stars 
form in densely packed clusters many form in widely spread OB associations. 

In any case, it is possible that a collisional build-up of high-mass stars can 
occur, especially for the most massive stars in very dense clusters. The problem 
is the very high stellar density of already massive or at least intermediate-mass 
stellar objects that is required to get the process going. We can estimate the 
threshold stellar number density Ugtar or, equivalently, the average star-star 
separation (s = n^^^^) in a dense cluster for stellar collisions (encounters) to be 
important. Because of the fundamental role of gravitational focusing coupled 
with a distribution of velocities (i.e., a finite probability of low relative velocities 
at infinity Vqo), the cross section Cgrav, where 



for two stars of mass passing each other at periastron within a minimum 
distance Rmin is vastly enhanced over the geometrical cross section. Gravi- 
tational focusing can enhance the effective cross section (by factors of 10^) , 
thus rendering close stellar encounters in very young stellar clusters realistic. Of 
course, the condition to meet for massive star growth through mergers is that 
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the collision (close encounter) time tcoii for the stars to collide must be shorter 
than the timescale for stellar evolution of the most massive star in the cluster 
(^3Myr). 

Using the formula for the stellar collision time per star (Binney & Tremaine 
1987, Dale & Davies 2006), i.e., 

1 



Tcoll 
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Here we have normalized the expression with reasonable numbers for the num- 
ber density, mass, size, and velocity dispersion of the stellar collision partners, 
assumed to be of equal mass. We ignore binary stars for the time being. Note 
that the velocity dispersion is the ID-velocity dispersion, which is 2.3 km s^^ in 
the Orion Nebula cluster (van Altena et al. 1988), but is ^5kms^^ in denser 
and more massive cluster cores such as NGC 3603 and R136, or even larger in 
young globular clusters such as those in the Antennae (Mengel et al. 2002). 

Equation 5 is based on equal mass collision partners; however, recently the 
formula has been extended to cover nonequal mass encounters (Moeckel & Bally 
2006, 2007), as well as larger cross sections due to circumstellar disks and binary 
components (Davies et al. 2006) . This can decrease the threshold stellar number 
density for collisions from 10^ or 10^ pc"'^ down to about lO^pc"'^ (Bonnell & 
Bate 2005). 

In Figure 13 we show 3D SPH numerical simulations of the collision process 
of two pairs of stars, one pair of equal mass and one pair of unequal mass. 
Interim Summary 

The primary difference between massive star formation by monolithic collapse 
versus by competitive accretion is that in monolithic collapse the mass is as- 
sumed to be gathered before the star-formation process begins, whereas in com- 
petitive accretion the mass is gathered during the star-formation process. If the 
former is true, we should find massive prestellar cores that live for long times. 
But how did they form in the first place? Perhaps turbulence could support 
a slow build-up, but turbulence induces density fluctuations, so why does the 
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quasi-static core not fragment into a cluster of low-mass stars? (This is the 
Dobbs, Bonnell & Clark 2005 versus Krumholz, Klein & McKee 2007 debate, 
see also Krumholz 2006). Magnetic fields may help to stabilize massive cores 
against subfragmentation. 

The strength of competitive accretion is that it provides a physical mech- 
anism to gather the matter. The gravitational potential of the protocluster 
clump or cluster of stars funnels a significant fraction of gaseous material to the 
cluster center, there to be accreted by the protomassive stars. In addition, the 
steep nonlinear mass dependence of the accretion rate onto a gravitating point 
mass - the Bondi-Hoyle rate - naturally gives a steep, power-law IMF, close 
to the observed IMF in young clusters (see Bonnell, Larson & Zinnecker 2007) . 
In this scenario, massive star formation is explained in the context of low-mass 
star formation (cospatial clusters). 

Stellar mergers will be rare and only relevant for the most massive stars in the 
richest young clusters (such as young globular clusters). One reason to invoke 
stellar mergers is that ongoing competitive accretion (mass loading) increases 
the stellar density in the cluster center (mass segregation at birth), potentially 
to the point where grazing collisions become unavoidable. Stellar collisions 
with small impact parameters might be the process that forms rapidly rotating 
massive stars and hence the progenitors of (slow, long-duration) gamma-ray 
bursts. 

4.4 Triggered OB Star Formation 

The classical model of triggered OB star formation goes back to Elmegreen & 
Lada (1977). The idea is that the ionization shock front of one group of massive 
stars provides the external pressure to compress adjacent molecular cloud layers, 
thereby inducing the formation of a new group of massive stars, which in turn, 
by the same process, induces the formation of another generation of massive 
stars and so on. The Elmegreen & Lada (1977) model was developed to explain 
the sequence of spatially distinct OB subgroups in nearby OB associations such 
as Orion OBI or Sco-OB2 (Blaauw 1964, 1991). It is important to note that 
the Elmegreen & Lada model did not predict the formation of the observed 
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coeval, low-mass T Tauri star population in the subgroups of OB associations 
(Preibisch & Zinnecker 2007, Briceno et al. 2007). The low-mass population 
was assumed to form independently and in many locations spread out over the 
cloud. Hence the Elmegreen & Lada (1977) process of sequential triggered star 
formation, contrary to accepted wisdom, may not be the main mechanism that 
accounts for the existence of OB subgroups. Rather, supernova triggering could 
be at work, assuming that a supernova blast wave can trigger both high-mass 
and low-mass stars at the same time (which would then explain the coevality of 
the high- mass and low-mass stars). Numerical simulations along these lines are 
just beginning (e.g., Melioli et al. 2006). 

However, there are several clear cases known in our Galaxy where expanding 
HII regions have swept up molecular gas at their periphery and in which new 
massive stars have formed or are about to form. The latter is indicated by 
luminous IR sources (e.g., Sh 104 and ROW 79, Zavagno et al. 2005; ROW 
108, Comeron, Schneider & Russeil 2005). This is indeed reminscent of the 
Elmegreen & Lada (1977) model of triggered sequential star formation (the so- 
called collect and collapse scenario, Elmegreen 1998; see also Whitworth et al. 
1994 and Dale, Bonnell & Whitworth 2007). However, according to Elmegreen 
& Lada (1977), it may be expected that in a two-stage starburst (cf. Walborn 
& Parker 1992) the second generation protocluster hosts only a small group of 
high-mass OB stars without the concomitant multitude of low-mass stars. Such 
a small N-body group, much like the Orion Trapezium system, would then be 
highly dynamically unstable and could help explain the occurrence of runaway 
OB field stars (Clarke & Pringle 1992, Allen, Poveda & Hernandez-Alcantara 
2004). An important implication of this scenario would be that the HH region 
of a big star cluster (e.g., NGC 3603) is not able to trigger the formation of a 
similarly massive second generation star cluster, but only a smaller mass cluster 
(e.g., IRS9 in NGC 3603; see Niirnberger 2003), potentially with a top-heavy 
stellar IMF; similarly, the R136 cluster in 30 Dor does not seem to trigger the 
formation of a new massive cluster but just a few small groups of embedded 
protostars (Rubio et al. 1998, Walborn et al. 1999, Brandner et al. 2001, 
Walborn, Maiz-Apellaniz & Barba 2002). 

A very interesting case is N81, one of the most compact HH regions in the 
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Small Magellanic Cloud (SMC) . N81 is isolated in an area of low extinction in 
Shapley's wing, 1.2 kpc away from the main body of the SMC. In contrast to 
other compact HII regions, typically located within or at the edge of giant HII 
regions, it appears that SMC N81 has been formed in isolation. The study of 
its stellar inventory shows that this 'high-excitation blob' is ionized by at least 
eight near-ZAMS O stars in an instantaneous burst (Heydari-Malayeri ct al. 
2002, 2003). The question here is, Where is the trigger? How did these massive 
stars form? 

Back to the Galaxy, another region of note is the Carina nebula powered 
by the most extreme grouping of massive stars in the southern Milky Way 
(Smith & Brooks 2007). Here ^ 65 O stars (including many in the clusters Trl4 
and Trl6, but excluding r] Car) provide a total of 10^^ erg of kinetic energy 
and a Lyman continuum luminosity of 10^^ photons s~^, creating a giant 
superbubble. Most of the bubble (seen as polycyclic aromatic hydrocarbons 
in 5piteer/MSX images) resides as atomic gas in the photodissociation regions 
and not in dense molecular clouds. The synchronized star formation around 
the periphery of Carina strongly suggests that star formation was triggered by 
stellar winds. The second-generation population appears to involve a cascade 
toward preferentially intermediate-mass and low-mass stars (Smith & Brooks 
2007), but the situation may change soon when r] Car and its siblings cxplod(^ as 
supernovae, rejuvenating massive star formation. The idea that it is supcrnovae 
rather than HII regions and stellar winds that stimulate wide-spread massive 
star formation is not new (see Herbst & Assousa 1977 for the Canis Majoris star 
formation region, see also Gerola & Seiden 1978 for galactic spiral arms, and the 
pioneering paper by Opik 1953). However, what is new is the evidence from the 
Upper Scorpius OB subgroup data that supernovae can trigger both high-mass 
and low-mass star formation in an OB subgroup at the same time (Preibisch & 
Zinnecker 2007). Supernovae need a critical distance from a molecular cloud - 
not too close and not too distant - to be an effective trigger for star formation 
because of momentum transfer of the blast wave onto prestellar cores (Vanhala 
& Cameron 1998, Vanhala et al. 1998). Other triggering mechanisms, like 
radiatively driven implosion of globules, also operate, but seem to be secondary 
processes, forming only small stellar groups rather than whole OB subgroups 
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with thousands of stars (for a review of a whole variety of triggering mechanisms, 
see Elmegreen 1998). 

The Carina star-forming region (diameter 150 pc, age 3 Myr, total IR-luminosity 
W'^Lq) may be the galactic analog of giant extragalactic HII regions such as 
NGC 604 in M33 (size 40arcsec or 140 pc), a very large OB association (see 
Maiz-Apellaniz, Perez & Mas-Hesse 2004). 

Closer to home, the Orion Nebula cluster is likely an example of triggered 
massive star formation. The appearance of the Orion A and B molecular clouds 
suggest an interaction caused by energy input from the Orion OB lab association 
subgroups (Blaauw 1991). Furthermore, the formation of the nearest massive 
star in the Orion BN/KL region, which hosts a very bright IR source, has prob- 
ably been triggered by the Orion Trapezium star's HII region pushing against 
the Orion Molecular Cloud right behind the Trapezium cluster (see Bally 2002). 

Recent Spitzer/IRAC Galactic plane observations (Churchwell et al. 2006) 
have revealed dozens of parsec-sized bubbles, formed by hot young stars in mas- 
sive star forming regions. Among the 80 or so ring- like structures produced 
by young OB stars, several show secondary bubbles on the rim of the primary 
bubbles, suggestive of triggered star formation. However, is this morphological 
evidence conchisivc, i.e., is triggering necessary to explain such spatial correla- 
tions? In an interesting paper, Dale, Clark & Bonnell (2007) attempt to address 
this and similar issues. They conduct SPH simulations to examine the difference 
between triggered and revealed star formation. They study the impact of irradi- 
ation by an external source of ionizing photons on a turbulent massive molecular 
cloud and compare the results (the number and type of stars formed) with a 
control simulation where the turbulent cloud evolves without the impact of the 
irradiation. They find that, although the external ionization has a dramatic 
effect on the morphology of the model cloud, its impact in terms of extra star 
formation is surprisingly minor; the feedback effects can be both positive and 
negative, accelerating the formation of some objects and delaying the formation 
of others. Only a few objects form that would otherwise not have formed, and 
the effect of induced star formation on the overall star formation efficiency is 
less than a factor of two. 
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4.5 Dynamical Evolution: Mass Segregation and Run- 
away OB Steirs 

4.5.1 Mass Segregation 

Massive stars arc often found near the centers of star clusters but not exclusively. 
Examples where massive stars are preferentially located in and near the cluster 
center include the Trapazium in the Orion Nebula cluster (Hillenbrand 1997) or 
the WR and O stars in the NGC 3603 cluster (Drissen et al. 1995). The issue 
is whether this is a birthmark (we call this effect prompt mass segregation) or 
whether this is an N-body evolutionary effect after birth (in which case, we call it 
dynamical mass segregation). The way to decide between these two possibilities 
is to estimate the dynamical time for massive stars to sink to the center of the 
gravitational potential from the half-mass radius of the cluster. If this timescale 
turns out to be too long, i.e., longer than the age of the cluster, dynamical mass 
segregation is ruled out and prompt mass segregation is indicated (e.g., Bonnell 
& Davies 1998, who infer that mass segregation in the Orion Nebula cluster 
must have been prompt and hence an important constraint for understanding 
massive star formation). 

4.5.2 Runaway OB Stars 

Runaway OB stars were defined by Blaauw (1961) as massive stars with radial 
velocities in excess of 40 km s^^ . Their high space velocities prompt the question 
of how they were accelerated to such kinetic energies. The original proposal 
involved an asymmetric supernova in a massive binary system that provided 
the kick to eject the companion from the system (Blaauw 1961). A strong 
prediction of this model is that the runaway massive stars should themselves be 
single. 

More recently, an alternative suggestion for the runaway phenomenon was 
put forward (Poveda, Ruiz & Allen 1967, Gies 1987). Here the idea is that 
massive stars in dense cluster cores (which is their preferential location) undergo 
dynamical three-body encounters, especially if a massive binary is involved. In 
such gravitational interactions, potential and kinetic energy are exchanged: The 
binary orbit shrinks, and the intruder extracts the potential energy and converts 
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it into the corresponding kinetic energy. It can also happen that an exchange 
reaction occurs so that one of the binary components (usually the lighter one, 
rather than the intruder) is ejected. Conservation of momentum leads to a recoil 
for each of the two system components, so they fly apart in roughly opposite 
directions. In this scenario, the runaway stars need not be single but can, in 
principle, be binaries, at least in 50% of the cases. 

Hoogerwerf, de Bruijne & de Zeeuw (2001) give evidence that both processes 
for runaway OB stars occur in nature: while ( Ophiuchi is a very high proper 
motion massive 09.5 star without an antipode in the opposite direction, and 
hence an example for the supernova scenario, another famous case, that of 
H Columbae and AE Aurigae, argues for the dynamical interaction scenario, 
where two hard binary systems likely had a close encounter in the 5Myr old 
NGC 1980 Orion cluster, in which two of the binary components were set off in 
opposite directions with speeds of around lOOkms"^. In their wake is t Ori, a 
tight, eccentric, massive binary (Gualandris, Portegies Zwart & Eggleton 2004), 
known to be classic, coUiding-wind X-ray binary. Other famous examples of 
massive runaway stars include A Cep, ^ Per, and a Cam (astronomy picture of 
the day on 24 Nov 2006). 

The dynamical ejection model for runaway stars would predict not only very 
fast ejection speeds, but also slower ones, as milder interaction events with 
wider, i.e., softer, binary pairs can occur (Kroupa 2000). This must be kept in 
mind when discussing the origin of massive field stars: with an escape speed of 
only 5-10kms~^, a massive O star with an age of around 3-5Myr could still 
travel 15-50pc from its cluster birth place. Some 40 nearby field O stars are 
known and at least half of them can be associated with a cluster origin (de Wit 
et al. 2004); for the rest the situation is unclear. It is interesting to recall that 
the spectroscopic and visual binary frequency among these runaway O stars is 
very low (Mason et al. 1998). 
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5 MASSIVE STAR FORMATION: 

BINARY AND MULTIPLE SYSTEMS 

Some of the most important clues toward understanding the formation of mas- 
sive stars, that have been neglected in previous reviews (e.g., Massey 2003), 
come from their high frequency of binary and multiple systems, together with 
an analysis of their properties (period distributions, mass ratios, orbital eccen- 
tricities); sec Section 2.4 and the discussion later in Section 6.2. Massive 
stars often come in hierarchical triples with an almost equal-mass close massive 
binary and a third more distant companion. One example includes ff^ Ori A 
in Orion's Bar (see Preibisch, Weigelt & Zinnecker 2001). Another is a Ori 
in the center of the namesake cluster, whose hierarchical configuration is even 
more complex (see Sanz-Forcada, Pranciosini & Pallavicini 2004). These close 
massive binaries likely play an important role in the evolution and age dating 
of starburst galaxies, as pointed out by van Bever & Vanbeveren (1998), and 
an equally important role in estimating supernova rates and the rates of high- 
and low-mass X-ray binaries (see Verbunt 1993; see also H.A. Kobulnicky, C.L. 
Fryer & D.C. Kiminki, submitted). Massive close binaries also bias and ob- 
fuscate measurements of the velocity dispersion in dense starburst clusters (see 
Bosch et al. 2001). 

It seems worthwhile and appropriate to summarize the various formation 
processes of massive binary and multiple systems. The following five processes 
could be relevant. 

5.1 Disk or Filament Fragmentation 

This is the most obvious mechanism; it is similar to suggestions of how to 
form low-mass T Tauri binaries from filament fragmentation (Zinnecker 1991, 
Bonnell & Bastien 1992, Monin et al. 2007). A filamentary geometry may 
play a key role in the fragmentation process, because the isothermal case is a 
critical one for the collapse of a cylinder: the collapse and fragmentation of a 
cylinder can continue freely as long as the temperature continues to decrease, 
but not if it begins to increase (Larson 2005). As for disk fragmentation, gas 
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thermal physics controls the non-linear outcome of gravitational instability in 
low- mass circumstellar disks (Durisen 2001). The fragmentation of massive 
disks has recently been studied analytically by Kratter & Matzner (2006). They 
found that these disks are unstable to fragment if they are cold enough, and 
catastrophically so when M > lO^'' Mq yr^-'^. This mechanism accounts only for 
initially wide binaries. Cold rotating Keplerian disks produce circular orbits; 
sub-Keplerian disk rotation or cold filaments tumbling end over end produce 
highly eccentric orbits. 6^ Ori C may be an example of the latter (Kraus et 
al. 2007). Note that fragmentation is different from fission, the splitting up of 
a hydrostatic, rapidly rotating body. Fission does not work for a compressible 
fluid (Tohline & Durisen 2001); it only leads to the ejection of spiral arms and 
torques that slow down the rapid stellar rotation. 

5.2 Accretion onto a Low-Mass, Wide Binciry 

This is a less obvious but potentially key mechanism to form tight massive 
binaries. Speculated upon by Maeder & Behrend (2002), it was worked out by 
Bonnell & Bate (2005) . These authors realized that the orbital separation of the 
components of close binaries is much smaller than the Jeans length at reasonable 
gas densities and temperatures. Therefore, they investigated the physical idea 
that the orbital separation of an initially wide binary would shrink by letting 
both components accrete and grow in mass; the final separation then depends 
on the orbital angular momentum of the accreted material, as follows: 

Let us consider the angular momentum A = M v x R of a binary system. 
Because of voc (M/R)^/^ (v is the Keplerian orbital speed), 

A OC M3/2r1/2 

where M is the instantaneous total mass of the binary system and R is the 
separation of the binary components. If the accreted material has zero net an- 
gular momentum, as would be expected if the infall were spherically symmetric, 
then A remains constant (A = const), implying that RocM~^, i.e., the binary 
separation should be a strong function of the binary mass. If, instead, the ac- 
creted material has constant specific angular momentum (A/M = const), then 
the total angular momentum will scale with the mass of the binary (AocM), 
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implying that RaM~^. Thus it is easy to see that accretion onto a binary sys- 
tem can significantly decrease its separation at the same time that it increases 
its mass (Bate 2000). In the numerical SPH simulations, the early evolution 
is well parameterized by the relation RcxM^^, intermediate between the above 
two cases, indicating that the binary is accreting some angular momentum with 
the mass but that the net specific angular momentum of the accreting gas is 
decreasing with time. The reason for this is that random velocities of gas at 
large radii tend to cancel each other out, resulting in accretion of a lower net 
angular momentum. With an RaM~^ relation, we infer that a low-mass (e.g., 
3M0) wide pair (e.g., separation of order 100 AU) will become a binary system 
with a semimajor axis of around 1 AU by the time the system has accreted up to 
30 Mq. That is, accretional growth in mass by a factor of 10 will make the sep- 
aration decrease by a factor of 100, or more if the growth in mass is more than 
a factor of ten. This may explain the observed very tight 0-star spectroscopic 
binaries (see Section 2.4). 

The component masses, even if unequal in the beginning, would tend to 
become equal later when higher angular momentum cloud material falls prefer- 
entially on the secondary component (Bate & BonncU 1997), as the lower mass 
secondary is, by definition, further away from the center of mass and would thus 
carry the higher angular momentum in the system initially. 

Another very interesting new model to explain the origin of close massive 
'twins' is proposed by Krumholz & Thompson (2006). They invoke mass trans- 
fer in close, rapidly accreting protobinaries; this always pushes the initial bina- 
ries toward mass ratio unity. Their model is superficially similar to the model 
of Bonnell & Bate (2005), but in fact the physical details are quite different, 
involving the swelling of protostars undergoing deuterium shell burning. 

5.3 Failed Mergers in Stellar Collisions 

If massive stars can form through stellar mergers, then near misses might some- 
times form tight and eccentric massive binaries. A necessary condition is that 
the kinetic energy of the quasi-parabolic encounter can be dissipated as tidal 
energy (Fabian, Pringle & Rees 1975) in a gravitationally focused grazing fly- 
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by (Zinnecker & Bate 2002, Dale & Davies 2006). Apart from the very high 
stellar number densities required for this process to be significant and frequent, 
the main problem is the extreme fine tuning of the collision impact parameter. 
Indeed, the cross section is a very small concentric annulus in impact parameter 
space that is much smaller than the central area inside the annulus relevant for 
massive mergers (F. Rasio, private communication). 

One way out is disk formation during the tidal disruption of a lower density 
star by a higher density star in a non-head-on collision; see the example of 
a grazing collision between a 3 Mq premain-sequence star and a 10 Mq main- 
sequence star in Figure 13. This can lead to subsequent disk-assisted capture 
of a companion, to be discussed extensively in the next subsection (the so- 
called 'shred and add' scenario, Davies et al. 2006). A variant of the failed 
merger scenario is the off-center collision of two protostellar cores with extended 
accretion envelopes (Stabler, Palla & Ho 2000). These protostellar envelopes 
offer a much higher interaction cross section and can provide the necessary 
orbital drag for collisional massive binary formation (cf. Silk 1978). 

5.4 Disk-Assisted Capture 

If massive stars form via disk accretion, then the large disk radii increase the 
interaction cross section considerably. This suggests that disk interactions with 
neighboring stars could assist in capturing binary companions (Bally & Zin- 
necker 2005). Although this mechanism has been found insufficient for solar 
mass stars with disks (Heller 1995, Boffin et al. 1998), recent SPH/N-body 
simulations by Moeckel & Bally (2006, 2007) convincingly showed that disk- 
assisted capture is much more efficient in a regime suited to massive stars (ca. 
20 Mq) with large disks (ca. 500 AU). We note that this process works particu- 
larly well in providing massive stars with lower mass companions in rather wide 
orbits (separation ca. 100 AU). The same authors also discuss the consequences 
of a mass-dependent velocity dispersion and of an initial mass segregation for 
the capture rates. Furthermore, they considered the long-term survival of these 
binaries in a dense cluster. 

The point is that massive binaries with different separations and mass ratios 
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form by different processes. Although accretion seems to be the only way to 
form the observed tight equal-mass binaries in young clusters, fragmentation 
and disk-assisted capture can form the wider, unequal-mass binaries. 

5.5 N-Body Dynamical Evolution 

Small-N groups of massive and intermediate-mass stars are seen forming in 3D 
SPH simulations of the collapse and fragmentation of gas-rich protoclusters (e.g., 
Bonnell, Bate & Vine 2003). Binary formation in these groups is common and 
occurs through dynamical three-body capture (different from tidal capture) . To 
begin with, a massive star typically has a lower mass wide companion. With 
time, an exchange interaction with a more massive third object occurs, thus 
forming an almost equal mass but still wide binary. The wide binary then 
shrinks ('hardens') as it takes up most of the binding energy of the small group 
when the lower mass group members get kicked and escape. The typical final 
outcome is a tight equal-mass massive binary, with a lower mass wide compan- 
ion (Bate, Bonnell & Bromm 2002). This scenario explains many aspects of 
massive binaries, including the average companion star fraction of 1.5 and the 
prevalence of short-period spectroscopic systems. Also, the timcscales for the 
N-body dynamics are short enough for the whole process to occur on or before 
the ZAMS (van Albada 1968, Aarseth 2003). 

However, it would seem that the frequency of massive close binaries should 
correlate with the overall stellar number density in young clusters, which is not 
what is observed. In fact, at face value the opposite - an anticorrelation - is 
observed (Mermilliod & Garcia 2001). We have no real explanation as to how 
this environmental effect fits into the N-body picture, unless most of the most 
massive binaries in the densest young clusters have managed to merge into a 
single object. Stellar mergers are thought to occur in old globular clusters (e.g.. 
Dale & Davies 2006), and this is one possible explanation for the presence of 
blue stragglers (younger stars with masses higher than the cluster turnoff mass) . 
Stellar mergers may therefore be even more important in dense, massive, young 
globular clusters, both now and in the past. 



Massive Star Formation 



61 



5.6 The Origin of Trapezium Systems 

Finally, we address the problem of how Trapezium systems of massive stars, 
such as the one in the center of the Orion Nebula Cluster {6^ Ori), likely came 
into being. In brief, the idea is as follows: numerical SPH simulations of super- 
sonic gravo-turbulent fragmentation of a protocluster cloud (1000 Mq) suggest 
that a collapsing cloud develops a few subclusters (star+gas systems), which 
subsequently merge into a single cluster entity. Each subcluster carries one 
most massive star (likely already part of a multiple). Hence the merging of 
subclusters will result in a central Trapezium-type system (sec Figure 14), as 
observed in the core of the Orion Nebula Cluster (sec Figure 15). Note that 
components Al and Bl of the Orion Trapezium arc spatially unresolved eclips- 
ing spectroscopic binaries; for a summary of the parameters of the Trapezium 
multiple stars we refer to Preibisch et al. (1999) and Schertl ct al. (2003). Note 
also that component B is itself a Trapezium-like system, indicating the hierar- 
chical nature of massive star formation. The dynamical evolution of Trapezium 
systems, including stellar ejections, is discussed by Allen, Poveda & Hernandez- 
Alcantara (2004), taking into account their multiplicity substructure. Future 
studies of Trapezium-type systems will likely concentrate on embedded systems, 
such as W3-IRS5 (Megeath, Wilson & Corbin 2005). 
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6 MASSIVE STAR FORMATION: DISCUSSION 

6.1 Disks and Outflows 

Disks and outflows are a general phenomenon in low-mass star formation that 
is explained by the accretion-ejection connection (Camenzind 1990, Ferreira 
& PcUcticr 1995). That is, disk accretion energy powers either a disk wind 
(Pudritz & Norman 1983, Pudritz ct al. 2007) or an X-wind (Shu ct al. 1994, 
1995). These in turn arc magnetically coUimatcd into a molecular jet. The jet 
then runs into the parent molecular cloud and local ISM, thereby accelerating 
entrained ambient molecular gas into a wider momentum-driven molecular flow. 
Textbook examples include the famous Hubble Space Telescope (HST) source 
HH30 (Burrows et al. 1996) and the most beautiful IR jets HH211 and HH212 
(McCaughrean, Rayner & Zinnecker 1994; Gueth & Guilloteau 1999; Zinnecker, 
McCaughrean & Rayner 1998; Lee et al. 2006). 

Similar but usually less collimated outflows have been observed around more 
luminous and, hence, more massive young stars (Shepherd & Churchwell 1996a, 
b; see the compilation by Beuther & Shepherd 2005). This has led to the notion 
that high-mass star formation is a scalcd-up version of low-mass star formation. 
However, as noted by Arce et al. (2007), no highly collimated outflow has been 
observed for high-mass star-forming regions exceeding 10"'' L© (corresponding 
to ^25M0). Likewise, circumstellar disks have been detected around young 
B-typc massive stars (Zhang 2005), but disks around the progenitors of 0-type 
massive stars in excess of 20 M0 have eluded detection (Cesaroni et al. 2007). 
This suggests that simple-minded claims of exact similarity between high-mass 
star formation and low-mass star formation are partly based on a casual use 
of the term 'high-mass star' and partly on wishful thinking. Yes, the rotating 
disks around IRAS 20216-^4104 (Cesaroni et al. 1997, 1999; Cesaroni et al. 
2005b), G192. 16-3.82 (Shepherd, Claussen & Kurtz 2002), Cep A (Patel et al. 
2005), or GL 490 (Schreyer et al. 2006) are wonderful examples of millimeter- 
interferometric observations. But the central sources are at best early B-type 
stars (below 20 M©) rather than massive O stars where radiation pressure and 
photoevaporation play a role in the formation process. The same seems to be 
true for the large silhouette disk in M17 (Chini et al. 2004), which is associated 
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with a molecular hydrogen jet (Niirnberger et al. 2007). 

Let there be no misunderstanding. It is likely that the majority of massive 
stars forms by disk accretion. But the issue is really whether stars of 50 100 Mq 
can also form by some other, more dramatic process like stellar or protostcUar 
collisions (BonncU, Bate & Zinnecker 1998; Stabler, Palla & Ho 2000). One 
could imagine that the reason why no disks and coUimated outflows have been 
found for the most massive stars in statu nasccndi is because their disks are 
being photoeroded at the same time the stellar mass grows by disk accretion. 
The combination of a powerful stellar wind and radiative acceleration clear the 
polar regions of material and magnetic fields, thus leading to poorer collimation 
for their molecular outflows. Maybe the outflow from AFGL 2591 (Figure 16) 
is such a case. 

The most messy outflow from a massive star-forming region is from the 
BN/KL region in Orion. The OMCl outflow (mass 10 Mq, velocity 30-100 km s'^) 
and the H2 flngers (Allen & Burton 1993, McCaughrean & Mac Low 1997, Kaifu 
et al. 2000) have been interpreted as the result of a powerful explosion that oc- 
curred in the center of OMCl within the last 1000 yr. It has been proposed 
that a dynamical interaction 500 yr ago, possibly leading to a merger, may have 
produced the OMCl outflow and H2 finger system (Bally & Zinnecker 2005). 
In this scenario, interactions with surrounding gas have decelerated the impul- 
sive outflow powered by the stellar collision by about a factor of two, thereby 
reconciling the 1000-yr dynamical timcscalc with the 500- yr timcscalc for the 
interaction. The detection of oppositely directed motions in radio sources IRc2-I 
and BN provides support for models in which the OMCl outflow was powered 
by a dynamical interaction between high-mass or intermediate-mass stars in an 
ultradense environment (Rodriguez et al. 2005a; see also Tan 2004). 

Another issue that needs to be appreciated in the context of massive outflows 
is the question: to what extent do the observed massive outflows trace the 
combination of stellar wind and intrinsic disk wind? The latter is able to entrain 
a lot of ambient gas, and care must be taken not to mistake this gas mass for the 
intrinsic wind mass. This is important because the accretion rate is estimated 
from the outflow rate, so if the outflow rate is overestimated, the accretion rate 
is overestimated too. 
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Somewhat analogously, one must also be careful when deriving the stellar 
mass from the himinosity of the source (e.g., using IRAS and submm data) when 
part of this himinosity is nonstellar and comes from accretion onto the disk and 
onto the (proto-)star (see Section 3.4.2 and Figure 5). Accretion occurs as 
long as there is outflow activity, and the highest accretion rate translates into 
the most violent outflow. Spatially unresolved binary and multiple systems 
compound the issue (see the next section). 

A final issue is the inferred nonlinear mass dependence of the disk accretion 
rate on the instantaneous central stellar mass (Norberg & Maeder 2000, Behrend 
& Maeder 2001). This is a consequence of a misunderstanding perpetuated in 
the literature. It dates back to the observed relation that the rate of mass 
outflow increases roughly linearly with the source luminosity (for a recent plot 
see Wu et al. 2005). If the outflow rate scales with the disk accretion rate and 
the luminosity varies with some power of the stellar mass (e.g., square), then we 
obtain M, ~ (Zhang 2005). However, it is a leap of faith to interpret this 
relation as one of instantaneous quantities rather than one regarding average 
quantities (the latter simply means that the formation of more massive stars 
requires higher average accretion rates, whereas the former would imply that 
the accretion rate keeps increasing with increasing central stellar mass). By 
completely ignoring the contribution of the accretion luminosity in this scenario, 
one implicitly assumes that a high luminosity implies a high core mass, as 
opposed to a high luminosity likely indicating that the accretion rate onto a 
lower mass (proto)star is high! 

A problem with an ever increasing accretion rate is not only how to suddenly 
stop the accretion, but also the observation that outflow activity is strongest 
during the earliest phases of star formation rather than during the final phases. 
In fact, the same problem exists for the turbulent core scenario of massive star 
formation (McKee & Tan 2003) where the accretion rate is predicted to increase 
finearly with time, implying that the most vigorous outfiows should occur toward 
the end of accretion rather than at the beginning. 
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6.2 Binciry Statistics and the Most Massive Bineiry Sys- 
tems 

The first spectroscopic binary survey among O stars (V < 7) was conducted by 
Garmany, Conti & Massey (1980). Tliey found tliat 24 out of their 67 O stars are 
spectroscopic binaries (36%). Some 15 are double-lined and 4 are single-lined, 
mostly with 1-10 day periods; the remaining 5 O stars are binary candidates. 
A later important work is the optical speckle survey of O stars carried out 
by Mason et al. (1998). This is a magnitude-limited sample (V<8) of 227 
bright O stars, observed at O.larcsec spatial resolution (15 new detections). 
In addition, the previously known spectroscopic binaries were included in a 
comprehensive analysis. We summarize their results as follows. Besides the 50 
known spectroscopic binary systems with periods less than 30 days and typical 
mass ratios of 0.5-1.0, there is a similar fraction of wide visual binaries (42 
systems with probable orbital motions, with periods larger than 100- 1000 yr, 
and typical mass ratios from 0.5 to 0.1). Clearly the distribution of mass ratios is 
different for wide binaries and for close binaries. However there is a huge gap in 
orbital periods between close and wide binaries. This is a selection effect: Such 
binaries have been hard to detect with past techniques. Interferometry will be 
able to find such binaries in the near future. Mason et al. (1998) estimate that 
the total O star (spectroscopic -|- visual) binary frequency is about 60%, but 
they suggest that when the above gap is filled, the frequency may well be 100 %. 
They also note that most binaries occur in clusters and associations and that 
binaries axe less common among field stars and especially among runaway stars. 

10 - 25 % of all O stars are runaway stars, whereas among the B stars this 
fraction is only about 2 % (Gies & Bolton 1986). The fact that the 0-star run- 
aways are 10 times more common than the B-star runaways implies that the 
secondaries of 0-star primaries (which are the ones that get ejected, Leonhard & 
Duncan 1990) should also preferentially be O stars. That is, many 0-0 binaries 
must have mass ratios skewed to unity (Clarke & Pringle 1992), in accordance 
with the observed high frequency of double-lined spectroscopic binaries among 
0-type stars in young clusters. However, the binary frequency for 0-type run- 
aways is generally lower than that in young clusters (Mason et al. 1998). A 
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possible explanation could be that the runaway recoil process produces a tight 
binary and a single massive star ejected in almost opposite directions, thus pre- 
dicting a spectroscopic binary frequency for runaway O stars that is half that 
of young 0-star clusters. 

Wolf-Rayet (WR) stars are evolved massive stars and, strictly speaking, do 
not belong to our main-sequence statistics. On the other hand, their progenitors 
were massive stars that lost much of their mass because of strong stellar winds. 
The mass loss rate and terminal wind speed depend on metallicity; more metal- 
rich stars lose more mass and lose it more quickly. The initial masses of WR 
stars cannot easily be determined from observations alone. Many new WR star 
binaries as well as mixed WR-OB star pairs have recently been discovered. The 
important question is whether the binary frequency of WR binaries and WR- 
OB star pairs depends on metallicity. We leave a discussion of WR-star binary 
statistics and related topics to Crowther (2007, in this volume). 

The 03 stars HD 93129A and HD 93205 in the Carina young star clusters 
Trumpler 14 and 16 were once thought to be the most massive single stars 
in the Galaxy (Taresch et al. 1997, Antokhina et al. 2000). In the case of 
HD 93129A this claim is based on the detailed non-local thermodynamic equi- 
librium analysis of the UV and hydrogen and helium optical spectrum of this 
0-type supergiant (Toff = 52,000 K). This analysis implies a very high bolomet- 
ric luminosity, 2.5 x 10^ Lq. The extreme dynamical stellar wind properties 
(M = -2OM0 Myr"-'^, v = 3200 km s^^) can be used to infer a stellar mass of 
13O±2OM0 (Taresch et al. 1997). However, Walborn (2003) and Nelan et al. 
(2004), using the Fine Guidance Sensor on HST, recently resolved HD 93129A 
into two components with a separation of 60 mas or 165 AU and a visual mag- 
nitude difference of 0.5 mag. This corresponds to component masses of about 
80 and 50 Mq. It is unclear how the new stellar parameters can be consistent 
with the high luminosity and the extreme wind properties, but this example 
shows strikingly how we can be misled in our conclusions if the binary nature 
of massive stars is unrecognized. 

As for HD 93205, this object consists of an 03V and an 08V noneclipsing 
pair with an orbital period of about 6 days. Because the inclination angle of the 
orbit is unknown, we only know that Mi • sin^ i = 29 M© and M2 • sin^ i = 13 M© 
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(Antokhina et al. 2000). An 08V star should correspond to about 2O-25M0, 
so one can tentatively infer a primary mass Mi = 50- 60 Mq. Morrell et al. 
(2001) suggest a spectral type 03.5 V and a mass of 48 Mq for the primary of 
HD 93205, making it one of the earliest main-sequence stars in the Galaxy at 
this time. There is at least one other massive double-lined spectroscopic binary 
in the Trumplcr 16 cluster (Tr 16-110 with P = 3.5 days), as well as two massive 
eclipsing binaries (Tr 16-1 and Tr 16-104) with P = 2.2 days (Rauw, personal 
communication). Also in the same cluster, there is one of the most massive 
WR stars in the Galaxy (WR 22, HD 92740), which is the primary of a spec- 
troscopic binary with P = 80 days. Its current mass is 45 Mq (Schweickhardt et 
al. 1999), corresponding to a ZAMS mass ~8OM0 (Rauw, personal communi- 
cation) . We refer to Walborn et al. (2002) for further discussion of the most 
massive stars known (02V stars). These include Gyg OB2-22A and Pismis 24-1, 
with masses in excess of 100 M© and, up until recently, believed to be single; 
however, the latter has now been resolved by HST as a visual triple system 
(see Mai'z-Apellaniz et al. 2007). A similar case was LBV 1806-20: Eikenberry 
et al. (2004) determined its mass to be close to 200 M©, but Figer, Najarro 
& Kudritzki (2004) dethroned it, finding that it is most likely a spectroscopic 
binary. Therefore, at present, the 'Pistol Star' with an estimated initial mass 
of 200- 250 M0 (Figer et al. 1998) in the Quintuplet Cluster near the Galactic 
Center is still considered to be the most massive star in the Galaxy. 

However, the most massive stars with dynamically measured masses reside 
in the double-lined spectroscopic, eclipsing binary WR20a: the primary and 
secondary component masses are 83 Mq and 82 Mq, respectively, with an error 
bar of 5 Mq . Rauw et al. (2004) originally discovered the spectroscopic binary 
nature of the object with an orbital period of about 3.7 days. Bonanos et al. 
(2004) obtained an eclipse light curve in the I-band and so could determine the 
inclination angle and hence the masses. Both objects are slightly evolved (SpT 
WN6ha), and still undergoing core hydrogen burning, but on their way from 
early 0-type to WR stars - with no previous phase of Roche lobe overflow. The 
binary system lies in the Westerlund 2 cluster associated with the RCW 49 HII 
region, but interestingly not near the cluster center. Speculation has it that 
this most heavy close binary system was ejected from the cluster core due to 
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a dynamical interaction. A spectacular 5piteer/IRAC image of the RCW 49 
HII region and cluster can be found in the Infrared Legacy Gallery at IPAC, 
courtesy of E. Churchwell and NASA/JPL-Caltech. A glimpse of some other 
spectacular Spiteer/IRAC images of high-mass star forming regions is given in 
Brandl et al. (2005), including the embedded proto-OB association W49A. 

Additional striking examples of multiplicity among very massive stars can be 
found in the 30 Dor cluster in the LMC, where Massey, Penny & Vukovich (2002) 
have identified four tight, double-lined spectroscopic binaries, three of which axe 
eclipsing systems (R136-38 is the most massive system, with component masses 
of 57 and 23 M©, while R136-42 is a close second with component masses of 
40 and 33 (the latter is actually a physical pair of 03V stars with a period 
of 2.89 days). Other massive stars in 30 Dor, including Melnick 34, are being 
monitored for radial velocity variations. 

Unlike the center of NGC 3603, which is dominated by a massive WR binary 
system, a recent, still unpublished near-IR integral field spectroscopic survey of 
the R136 central regions with the ESO-VLT has not identified any WR spectro- 
scopic binaries (T. Mofi^at & O. Schnurr, personal communication). This result 
is contrary to the physical intuition that would lead us to expect the heaviest 
objects, i.e., massive binaries, to be located in the very center (mass segregation 
due to dynamical friction). 

The center of the R136 cluster that was once considered to host a super- 
massive star (Feitzinger et al. 1980; Cassinelli, Mathis & Savage 1981) actually 
consists of a group of eight massive stars within a projected radius of 0.4arcsec 
or equivalently 0.1 pc (Weigelt & Baier 1985; Pehlemann, Hofmann & Weigelt 
1992), corresponding to a central mass density in massive stars of 10^ Mq per 
cubic parsec. We do not know how much mass there is in low-mass stars in the 
center, as two-body dynamical mass segregation could have led to a depletion 
of low-mass stars there. The density of massive stars alone may not be high 
enough at this time to enable a runaway stellar collision process, leading to the 
formation of an intermediate mass black hole, as envisaged by Portegies Zwart 
& McMillan (2002). It is enough, though, for the occasional collision between 
two massive stars, leading to the formation of a rapidly rotating massive object, 
possibly the progenitor for a gamma-ray burster (Fryer & Heger 2005, Zinnecker 
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2006b). Although the R136 cluster is quite dense in its center, more massive 
young globular clusters such as those seen in the Antennae merging galaxies 
may be denser still, and a runaway collision process cannot be excluded. 

6.3 The Universality of the Upper Initial Mass Function 

How can we understand the universality of the upper IMF? Wherever we look, 
the data seem to be compatible with a Salpeter power-law with a logarithmic 
slope of -1.35. Even population studies at moderate to high redshift are consis- 
tent with the assumption of a Salpeter power law (e.g., Baldry & Glazebrook 
2003). 

In essence, two completely different schools of thought have attempted to ex- 
plain the robustness of the Salpeter IMF despite various kinds of environmental 
factors that could be expected to change this distribution (such as metallicity, 
gas pressure, or the density of stellar systems, to name but a few). These are 
competitive accretion and random sampling of fractal clouds. 

As for competitive accretion, this has been discussed extensively in the re- 
view of Bonnell, Larson & Zinncckcr (2007) and will not be repeated here. 
Suffice it to say that competitive accretion of protostars for cluster gas can ex- 
plain the observed Salpeter stellar mass distribution for massive stars, if most 
massive stars form in dense clusters (Bonnell ct al. 2001b, Klcsscn 2001a, b, 
Klessen & Burkert 2001; for a discussion of the onset and richness of clustering 
as a function of the most massive cluster star see Testi et al. 1997 and Testi, 
Palla & Natta 1999). 

If, however, a substantial fraction of massive stars forms in less dense OB 
associations (Garmany 1994; Clark, Bonnell & Zinnecker 2005 ) or in isolated 
places in the field (for example, where was Betelgeuse - a ~ 20 Mq red supergiant 
- born?), then competitive accretion cannot be invoked, and another idea is 
needed (e.g. Lamers et al. 2002; Li, Klessen & Mac Low 2003). 

Elmegreen (1997) presented such an idea, essentially a geometrical model, 
with a physical icing on the (geometrical) cake. This model is based on a random 
sampling of mass in turbulent fractal interstellar clouds (see Figure 17). A 
star was assumed to get a fixed fraction of the gas mass of the cloud piece in 
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which it formed. The mass distribution of the pieces in any hierarchy has a 
logarithmic slope of -1 (i.e., equal mass in equal logarithmic intervals), so the 
IMF would seem to end up with this slope, but the sampling rate for pieces 
was assumed to be proportional to the square root of the local gas density to 
mimic the dynamical processes that are involved. This local density depends 
on the level in the hierarchy according to the fractal scaling of density with size 
(because mass scales with size to a power equal to the fractal dimension, density 
scales with size to a power equal to the fractal dimension minus three). This 
density dependence means that lower mass regions are sampled more frequently, 
and the mass function slope steepens from -1 to -1.35, which is the Salpeter 
slope. The net result of this sampling is a mass function for model stars that is 
indistinguishable from the observed IMF for young clusters and OB associations. 

Another interesting attempt to connect star formation to fractal cloud struc- 
ture was made by Henriksen (1991, see his figure 1), but the details differ from 
Elmegreen's (1997) model. An important variant of the fractal molecular cloud 
IMF model is the concept in which the stellar mass function derives directly 
from the mass distribution of cloud cores. Recent observations show a surpris- 
ing similarity between the scaled cloud core mass distribution and the stellar 
IMF (Alves, Lombardi & Lada 2007). 

Yet another alternative for explaining a power-law IMF is due to Basu & 
Jones (2004). They note that the power-law tail in the mass function of pro- 
tostellar condensations and stars arises from the accretion of ambient cloud 
material on to the condensation, coupled with a nonuniform (exponential) dis- 
tribution of accretion lifetimes (cf. radioactive decay). Thus, this model assumes 
that not all condensations accrete for the same time. If we start with proto- 
stars with a log-normal mass distribution (possibly a reasonable assumption 
expected from the central limit theorem, see Zinnecker 1985), this log- normal 
distribution develops a power-law tail at high masses if the accretion rate is 
directly proportional to the instantaneous mass of the accreting object and if 
the probability of stopping accretion is constant in time. The latter implies 
an exponential probability distribution of accretion timescales with a constant 
death rate. How general this random accretion model is and how it relates to 
the competitive accretion model above needs to be further explored (cf. Bate 



Massive Star Formation 



71 



& Bonnell 2005). 

Needless to say that the characteristic mass scale of fragmentation and the 
IMF (a few tenths of a solar mass) is beyond the scope of this review, and we 
refer to Larson (1985, 2005) and also Whitworth, Boffin & Francis (1998) for an 
in-depth discussion of this issue. Finally, we refer to Klessen, Spaans & Jappsen 
(2007) for the first numerical hydrodynamical calculations of the characteristic 
stellar mass in starburst regions, predicting a top-heavy IMF. This is naturally 
explained as a consequence of the elevated thermal Jeans mass in the warmer 
and dustier starburst environment (including the Galactic Center). 

6.4 The Number of Accreting Massive Protostars in the 
Galaxy 

In this subsection, we derive an estimate of high- mass accreting stellar objects in 
the Galaxy (note that we arc trying to avoid the term protostars). To this end, 
we first calculate how much stellar mass is formed in the Galaxy over a time- 
interval tacc corresponding to the accretion phase (200,000 yr). Using a total 
Galactic star-formation rate (SFR) of 5M0yr~^ (Smith, Mezger & Biermann 
1978, Diehl et al. 2006) and multiplying by 200,000 yr, we obtain IO'^Mq. Of 
these, we assume a fraction fw 0.10 ends up in massive stars (above lOM©, say). 
The rest goes into intermediate- and low-mass stars according to a Salpeter 
(1955) field star IMF, with an effective lower mass limit of 0.1 Mq and upper 
mass limit of 100 M© (see Smith, Mezger & Biermann 1978, their table 2). 
A more realistic low-mass field star IMF (Kroupa 2002, Chabrier 2003) would 
increase the fraction f by about a factor of 2, yet would decrease the SFR derived 
above by about the same factor, thus leaving the product f x SFR approximately 
constant. 

The number (>M*) of stellar progenitors in the Galaxy to be found in 
the accretion phase, which eventually attain a final mass above a given value 
M*, can then be derived from the two equations (the first solving for Nq), 
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No / m-'^dm = f x SFR x Atacc 

N*(> m*) = No / m'^'dlogm 

Here m» =M»/M0 is the dimensionless stellar mass and dN/dlog(m) =Nom~'^ 
is the power-law upper stellar IMF with a slope x (x= 1.35 for a Salpeter IMF). 
Solving this equation results in the numbers N* (> M*) given in Table 4. 

Table 4: Expected number N of stellar progenitors in the Galaxy to be found in 
the accretion phase", which eventually attain a final mass above M*, for various 
logarithmic slopes x of the IMF** 



M* 

[M©] x=1.3 x=1.35 x=1.7 



> 10 


6300 


5400 


1600 


> 20 


2400 


2000 


480 


> 30 


1300 


1000 


220 


> 50 


480 


390 


74 



"Assumed to last 200,000 yr. 

''Assuming a total Galactic star formation rate of SM© yr~^. 

The table shows that we can expect about 1000 accreting massive stars in the 
Galaxy with masses in excess of 30 M© (i.e., spectral types earlier than 06.5V, 
see Table 3). Assuming lOkpc as the reference radius of the Milky Way, we see 
that the average surface density of accreting early 0-type massive stars is about 
3objects/kpc^. This demonstrates the rareness of massive protostars, but is in 
agreement with the few luminous and totally embedded objects known, such 
as BN/KL and Cep-A (see Table 5 in the Appendix which provides a list of 
massive star forming regions within 1 kpc from the Sun and some properties of 
the exciting stars). 
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6.5 Is There a Mciximum Stellar Mass Set by Steir Forma- 
tion? 

In Section 2.5 we reviewed the observational evidence for an upper stellar mass 
limit and concluded there was such a limiting mass somewhere in the range 
between 100 and 200 Mq , at least in Population I systems. Here we present a 
few thoughts on the possible physics of the stellar upper mass limit. 

The first question is whether this limit is set by stellar stability considera- 
tions or star-formation theory. Ledoux (1941), and later Schwarzschild & Harm 
(1959), concluded from linear stability analysis of radial stellar pulsations that 
there is an upper mass limit of the order of 100 Mq, corresponding to the ex- 
istence of a vibrational instability owing to nuclear reactions. This refers to 
massive stars of Population I chemical composition, whereas for lower chemi- 
cal abundances the unstable mass is higher. Appenzeller (1970) and Ziebarth 
(1970) concluded, from nonlinear calculations, that this vibrational instability 
(also known as e-instability) does not have a global disrupting effect, except 
perhaps in stars with masses greatly exceeding the critical mass. A more recent 
discussion can be found in Stothcrs & Chin (1993). These authors concluded 
that, for a mctallicity equal to or greater than 0.02, no nuclear-induced pulsa- 
tional instability developed for masses up to at least 150 Mq. For a metallicity 
equal to 0.004, corresponding to the value of the Small Magellanic Cloud, the 
critical mass for such an instability was still very high, around 140 M©. 

Larson & Starrfield (1971) were the first to stress a possible upper mass limit 
resulting from star-formation theory. They came up with a value of 50 - 100 M© 
as a larger and larger fraction of the growing mass is thrown out by radiation 
before the star reaches stellar conditions. 

Today, we realize that the opacity of dusty gas does not determine the upper 
mass limit, whereas the opacity of ionized gas (electron scattering and UV 
line blanketing) certainly may. The constraint imposed by dust opacity (see 
Figure 6) shows that a luminous object with an L/M ratio consistent with 
a main-sequence star M* > 100 M© can still accrete material as long as the 
object appears to be cool (<50K). This allows material to fall onto the disk. 
Within the disk much of the radial radiation flux is diverted into the polar 
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direction (flashlight effect). Provided angular momentum transfer is adequate, 

dusty material can thus flow radially inward until at some point the dust is 
destroyed (sublimation radius). Dust-free neutral gas has a lower opacity and 
it can continue to flow inward. Eventually, as the disk gas gets closer to the 
accreting star, the gas is ionized and the opacity increases again sharply. Disk 
models currently do not allow an accurate estimate of the opacity of the ionized 
disk gas close to the star; the lower limit to the opacity provided by electron 
scattering limits the maximum L/M to ~ 6 x 10^, well in excess of that expected 
from a main-sequence star of M, = 200 M©. 

Because very massive stars also possess strong radiation-driven winds, it 
is clear that the opacity of ionized stellar gas is sufficiently high that radiative 
acceleration exceeds surface gravity. The upper mass limit could thus result from 
the fact that mass loss from the star matches or exceeds accretion. As the line- 
driven wind mass loss from massive stars is metallicity-dependent (Kudritzki 
2002) the upper mass limit would be metallicity-dependent, too. However, if 
mass loss is by a continuum-driven instability, as suggested for jj Car and other 
Luminous Blue Variables (Smith & Owocki 2006), the upper mass limit would 
be independent of metallicity. 

It is unclear if stellar collisions and mergers in dense stellar systems can 
beat the opacity-limited accretion discussed above and form stars with masses 
beyond 200 M©. The merging of a binary star with two lOO-M© components 
is conceivable (Zinnecker 1986, Bonnell & Bate 2005), doubling the single star 
upper mass limit. It is likely that a coalescence process would deplete the 
number density of stars to collide with near the center of a deep cluster potential, 
thus imposing a density-dependent limiting mass, especially for young clusters 
undergoing core collapse (Portegies Zwart et al. 2006). Interestingly, in this 
case, too, the upper mass limit would be independent of the heavy element 
abundance. We refer to Omukai & Palla (2003) for a theoretical prediction of 
the increased upper mass limit (~6OOM0) in zero-metallicity. Population III 
conditions, albeit based on ID models. 
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6.6 Evolutionary Sequence 

We can characterize the formation of massive stars by the following crude four- 
stage evolutionary sequence: CDMC — > HDMC — > DAMS — > FIMS. Here 
CDMC means cold dense massive core, HDMC means hot dense massive core, 
DAMS stands for disk-accreting main-sequence star, and FIMS denotes the 
final main-sequence star. The FIMS star is more evolved than the theoretical 
concept of the ZAMS star, which assumes no prior hydrogen consumption. Of 
course, the sequence needs a more detailed description. The CDMCs are starless 
turbulent gravitationally bound condensations; they have near-virial equilibrium 
and are either on the verge of collapsing or, indeed, already collapsing. In the 
center an intermediate-mass protostar forms that can heat up the dense massive 
core, turning it into a HDMC. At this stage, coUimated jets and outflows flrst 
appear, traced by H2O and later by methanol maser emission. The central star 
grows in mass primarily by disk accretion and quickly becomes a DAMS object 
powered more by hydrogen burning than by disk accretion. The coUimated 
outflows become less coUimated with widening opening angles. At this stage, 
the accretion disk starts to get photoionized and partly photoevaporated, giving 
rise to a gravitationally confined hypercompact HII region (HCHII) with broad 
hydrogen recombination lines. 

Eventually the star has accreted its final mass and most of the disk has 
been dissipated. At this point the ionizing radiation is no longer quenched by 
the accretion flow and can expand freely into the vicinity of the star, blending 
together with ionized bubbles from other similarly massive stars. This is when 
we speak of an ultracompact HII region (UCHII), often formed by a group of 
ionizing stars within a volume of (0.1 pc)^ (e.g., W3-IRS5). UCHII regions 
frequently exhibit OH maser 18-cm radio emission [e.g., W3(0H)], which can 
be used to obtain a magnetic field strength by measuring the Zeeman splitting 
of the 1667 MHz line. Subsequently, the UCHII region keeps expanding and 
evolves into a compact radio HII region and finally into a normal diffuse optical 
HII region like the Orion Nebula. 

How this simple-minded picture changes when we deal with a massive star- 
burst protocluster, like the progenitor of R136 with dozens of massive stars 
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born in a small volume, is anybody's guess. It is conceivable that competitive 
accretion takes place in a gravitationally bound dense HII region. It is even con- 
ceivable that competitive accretion turns into cooperative accretion (E. Keto, 
private communication) which means the ionized flow pulled in by the collec- 
tive gravity of the cluster is redirected onto the less massive members of the 
cluster. Instead of being greedy and competing for additional mass with their 
siblings, the near-Eddington-limit massive stars start to deflect the inflow onto 
neighboring stars. By doing so, they aid the further growth of their lower mass 
neighbors. Thus, accretion becomes cooperative! Another possibility is the hi- 
erarchical merging of subclusters which naturally leads to prompt initial mass 
segregation in the resulting final cluster (McMillan, Vesperini & Portegies Zwart 
2007). In this scenario, many but not all subclusters (gas-|-stars) merge with 
each other and sink to the center of gravity owing to the loss of kinetic energy 
in these highly inelastic collisions. There they form a dense stellar cluster of 
massive stars, which will expel the residual cluster gas. Some subclusters miss 
and won't collide; these will form massive stars in the cluster outskirts. Such 
a core-halo structure is observed in many young clusters, including R136 (e.g., 
Moffat, Drissen & Shara 1994). Super OB associations like NGC 604 in M33 
have more widely spread subclusters that are unable to merge, thus avoiding 
dense cluster formation (cf. Hunter 1995). 
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7 WHY HIGH-MASS STAR FORMATION IS 
NOT A SCALED-UP VERSION OF LOW- 
MASS STAR FORMATION 

In this section, we try to convince the reader that high-mass star formation 
is indeed different from low-mass star formation, and not merely a scaled-up 
version of star formation by disk accretion for higher accretion rates. Many 
complex new physical processes enter the scene for high-mass stars. 

To start with, radiative forces on gas and dust play little or no role in the 
build-up of low-mass, solar-type stars, whereas a substantial fraction of the lu- 
minosity of high-mass stars is emitted in ionizing radiation, which introduces 
new effects such as the photoevaporation of the star's accretion disk and pro- 
tostellar envelope. This dramatically limits late accretion and the final stellar 
mass. In addition, the ionizing photons can photoevaporate the disks of the 
neighboring lower mass stars (cf. the proplyds in the Orion Nebula). The 
nonionizing far-UV radiation will influence the massive stars' molecular cloud 
environment by dissociating H2 and CO molecules. This requires photons of 
about 11.2 cV; even early-type B stars can produce these photons, but low- mass 
and intermediate-mass stars cannot. The radiative acceleration of dusty and 
gaseous matter also leads to radiation-driven bipolar winds, and ionizing radi- 
ation can escape through these wind-blown cavities (flashlight effect); however, 
the bipolar outflows from low-mass stars are generated by magneto-centrifugal 
forces. 

The second big difference between the formation of low-mass and high-mass 
stars is the fact that massive stars are practically born on the main sequence, 
whereas low-mass stars spend a considerable part of their youth as contract- 
ing premain-sequence objects (30Myr for a solar-mass star). A massive star, 
forming by accretional growth from an initially low-mass star with an accre- 
tion rate of lO~*M0yr~^ begins central hydrogen burning after about 9Mq 
have accumulated - 13 Mq for an accretion rate of 10~^ Mq yr~^. For low-mass 
stars, circumstellar disk evolution proceeds during the whole extended premain- 
sequence phase, whereas for massive stars the disk lifetime is very short (less 
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than IMyr). When we see massive O stars close to the ZAMS (Walborn 2007), 
their disks have been dissipated (but see Kastner et al. 2006). It is only in the 
embedded phase that circumstellar CO-bandhcad emission and hydrogen emis- 
sion lines are detected (Bik ct al. 2005, Blum 2005), indicating a dense neutral 
remnant disk with an ionized upper layer. Bipolar outflows are another strong 
indication of the existence of disks. 

A further significant difference between low- and high-mass star formation is 
the role of competitive accretion in protoclusters. Competitive accretion is far 
more important for high-mass stars than for low-mass stars. The latter can form 
directly by Jeans-type gravitational instability and turbulence-induced cloud 
fragmentation. The former must accrete large amounts of protocluster gas. 

Gravitational dynamics (N-body interactions) also have a much greater effect 
when massive stars are involved. This can best be seen when considering the 
dynamical ejection of members of massive multiple systems producing runaway 
OB stars. This phenomenon is largely absent for low-mass stars. It appears 
that massive stars, when they form in their own local subclusters, are always 
accompanied by a small group of lower mass stars, some of which stay bound to 
their big parent star even after subclusters merge. Thus the higher companion 
star fraction observed for massive stars compared to solar-type stars testifies 
to a more dynamic scenario for the heavy objects thanks to their higher than 
average gravitational attraction. 

Finally, massive stars have a much bigger influence on triggering new star 
formation in adjacent regions than low-mass stars. Massive stars provide exter- 
nal pressure in the form of expanding HII regions, stellar winds, and supernovae 
explosions. They are capable of sustaining sequential and self-propagating star 
formation (Gerola & Seiden 1978), a process that low-mass stars are incapable 
of. Through runaway OB stars, massive star formation can trigger further mas- 
sive star formation over large (kpc) distances, an important feature in sustaining 
large-scale nuclear starbursts. We speculate that, if massive star formation is 
massively triggered, the individual collapse of massive cores is outside-in, instead 
of inside-out (Banerjee, Pudritz & Anderson 2006), and always magnetically su- 
percritical (Shu, Adams & Lizano 1987). In other words, magnetic fields likely 
play a more passive role in massive star formation, whereas in low-mass stars 
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this is opposite. 

In defense of the notion of a scaled-up formation picture, some arguments 
have been raised. Foremost is the observation that coUimated outflows occur 
both in high-mass and low-mass young stars (Beuther et al. 2002, Davis et al. 
2004). However, a closer look reveals that coUimated outflows and jets do not 
occur in the most luminous sources, but only in sources with total luminosities 
up to 10'' Lq (Shepherd 2005), with one or two exceptions (Garay et al. 2003, 
Rodriguez et al. 2005b). This implies that massive star formation can be a 
scaled-up version of low-mass star formation, but only up to early B stars. The 
outflow morphology does not scale further for O stars, which generate powerful, 
wide-angle, ionized winds, calling into question the relationship between outflow 
and accretion (Shepherd 2005). 

Another point that has been made relates to the IMF. It has been suggested 
(see the discussion in Zinnecker 2004) that there should be a feature (knee) in 
the IMF, at a critical mass where the intermediate-mass and high-mass star- 
formation processes diverge. However, this feature has not been seen and the 
upper IMF is a power-law with a constant slope (Kroupa 2002). It is unclear how 
to refute this argument, but the cxHistant slope of the upper IMF may simply 
imply that stellar collisions and mergers are not a dominant factor in massive 
star formation, except for the most massive stars (Bonnell & Bate 2002) where 
statistical fluctuations can hide a slope change. CoUisional growth of massive 
stars (Bonnell, Bate & Zinnecker 1998; Bally & Zinnecker 2005) also seems to 
be ruled out by recent observations of stellar rotation: There is a continuous 
behavior of the specific stellar spin angular momentum over the full range of 
stellar masses, i.e., J/M M*^-^ between 0.2 Mq and 50 Mq, indicative of a single 
stellar formation and angular momentum regulation mechanism (Wolff et al. 
2006). However, a coUisional process and stellar mergers are not ruled out for 
the most massive early 0-type stars (M > 50 M©). 
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8 OUTLOOK: RELEVANT FUTURE OBSER- 
VATIONS 

Here we suggest a number of future observations that would help promote?advaiice 
our understanding of massive star formation. The list, however, is by no means 
complete. 

1. IR-observations with an Extremely Large Telescope (ELT) 
(diameter 30 -42 m) 

As described in Zinnecker (2006a), such a powerful telescope could pene- 
trate the dust extinction of ultracompact HII regions (Ay = 100 - 200 mag) 
in the near-IR (Ak = 10 - 20) and see the stellar photospheres of mas- 
sive stars, resolving very tight embedded clusters, such as W3-IRS5 (cf. 
Megeath, Wilson & Corbin 2005) . There is hope that an ELT might have 
sharp enough vision to test the prediction (Zinnecker 2006b) that gamma- 
ray bursts occur in the centers of young massive protoglobular clusters. 

2. Submillimeter observations with the Atacama Large Millimeter Array 
The spatial resolution of the Atacama Large Millimeter Array (ALMA) 
at 350 microns (^ 0.01 arcsec, i.e., 50 AU at 5kpc) would allow us to map 
dense molecular clumps, where massive stars are believed to form, in the 
dust continuum and see whether fragmentation is occurring (Dobbs, Bon- 
nell & Clark 2005) or not (Krumholz 2006). ALMA should also be able 
to measure rotation curves of circumstellar disk gas around massive stars 
and thus measure their enclosed stellar mass (M. Krumholz, R.L Klein & 
CF. McKee, submitted). 

3. Mid-infrared observations with the James Webb Space Telescope 

The James Webb Space Telescope will have unparalleled background-limited 
sensitivity to reveal very embedded populations around young massive 
stars in the making in near- and far-galactic star forming regions as well 
as regions like 30 Dor/R136 in the Magellanic Clouds. Progress may also 
be expected for understanding star formation in the Galactic Center Re- 
gions, including resolution of the Arches and Quintuplet clusters and the 
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OB populations near Sgr A* (Figer 2003, Kim et al. 2006, Paumard et al. 
2006, Martins et al. 2006). 

4. Proper motion observations with Gaia 

Precise proper motion data of young clusters and OB associations are ea- 
gerly awaited, surpassing the limited milli-arcsec precision of Hipparcos 
(de Zeeuw et al. 1999). This will be crucial to understanding the dy- 
namical nature of OB associations, such as the question of whether an 
OB association is just an expanding, dissolving massive star cluster or 
something of its own (Blaauw 1983, Bricefio et al. 2007). Hence Gaia 
will tell us if the stars of an association expand from a common center or 
from several centers or from no center at all. Gaia will also provide us 
with a complete sample of field O stars and their parallaxes within a few 
kiloparsecs. 

5. Far-IR observations with Herschel 

The immediate future will see the launch of ESA's Herschel space ob- 
servatory, whose Photodetector Array Camera & Spectrometer (PACS) 
instrument will make it possible to follow up on Spitzer/MlPS observa- 
tions and to detect highly embedded massive protostellar sources down 
to a spatial resolution of 3 arcsec at 60 microns. Surveys with the Spec- 
tral and Photometric Imaging Receiver (SPIRE) will allow us to identify 
and characterize pre-coUapse cores at an early evolutionary stage (Boss 
& Yorkc 1995). PACS together with SPIRE can be used to better con- 
strain the 'collect and collapse' model of triggered massive star formation 
(Zavagno et al. 2005), discriminating between true triggering versus the 
compression of pre-existing clumps. 

6. Interferometric near-IR observations with VLTI and LBTI 

Very high angular resolution observations with closure phases will give us 
interferometric images (and not just visibilities) of ultracompact HII re- 
gions and the stellar content of hot molecular cores. This will allow us to 
obtain a clearer picture of the maternities of massive stars (Garay 2005), 
including images of disks and jets as well as astrometric information, such 
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as the stellar velocity dispersion, dynamics of the stars in protocluster cen- 
ters, and the orbital motion of short-period embedded binary and multiple 
systems. 

7. Long-term variability surveys The formation and evolution of massive 
stars is likely to involve a variety of cataclysmic or eruptive phenom- 
ena on short timescales (days to years), such as episodic FU Orionis-like 
accretion (Chini, private communication) and stellar mergers (Bally & 
Zinnecker 2005), or i] Car- like LBV eruptions (Smith & Owocki 2006) and 
even infrared supcrnovae in luminous starburst galaxies (Maiolino et al. 
2002). Without continuous monitoring, these (embedded) events may go 
unnoticed. This calls for long-term variability compaigns, in the near- 
infrared, for example from a medium-sized survey telescope at Dome C in 
Antarctica. 

8. "Observations" with a Numerical Telescope 

As numerical simulations of the star formation process become more re- 
fined in 3D resolution and the microphysics involved (magnetic fields, 
chemistry, ionization, heating, cooling, dust physics, etc.) the ability to 
accurately "observe" the evolution at various stages with respect to dust 
continuum features, line intensities and line profiles increases in impor- 
tance. In other words, theoretical modeling will have to keep pace with 
the rapid development of observing capabilities as described above. 
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9 FINAL SUMMARY 

Massive star formation is not a simple scaled- up version of low-mass star forma- 
tion, particularly when it comes to ZAMS O stars (>2OM0). The formation 
of early- type B stars, however, may still be considered a continuation of low- 
mass star formation. In terms of formation processes, the monolithic collapse 
and disk accretion model and the competitive accretion scenario axe the two 
opposite ends of a continuum of cases in the accretion theory of 0-type massive 
star formation. The coalescence process, introduced to circumvent the obsta- 
cle of radiation pressure, is no longer generally necessary but may still occur 
in exceptional circumstances, especially for very massive stars in the centers 
of very dense protoclusters or subclusters with strong initial mass segregation. 
The observed multiplicity and clustering of massive stars suggests complex N- 
body and gas dynamical interactions (tidal, drag, and capture effects) among 
the youngest stars or protostars that would contribute to the argument against 
a monolithic collapse of isolated massive protostellar cores. 

Massive stars are rare and correspond to the tail of a power-law stellar mass 
distribution. For a Salpeter IMF, the number of all newborn ZAMS OB stars 
(>8Mq) is a mere 10% of all stars in the mass range of 1 - 2Mq. At present, 
competitive accretion (i.e., runaway growth of a few objects from cluster gas) 
provides the best physical explanation of the high-mass IMF and its surprisingly 
universal Salpeter slope. The maximum stellar mass (around 150 Mq), long 
believed to be a result of stellar instability (Eddington-Ledoux limit) , might be 
due to the negative feedback caused by the increasingly destructive radiative 
erosion of massive accretion disks at increasingly higher stellar masses. This 
self-limiting star formation can, however, be beaten by stellar collisions under 
very special circumstances. 

Rapid external shock compression (i.e., supersonic gas motions) generating 
high column densities in less than a local free-fall time rather than slow quasi- 
static build-up of massive cores may be the recipe to set up the initial conditions 
for local and global bursts of massive star formation. Such bursts can both en- 
hance or quench further massive star formation (positive and negative feedback), 
depending on environment (gas density and pressure), but such a discussion is 
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beyond the scope of this review. The question of how feedback from massive 

stars can influence proto galaxy evolution and morphological types (cf. Sandage 
1986, Elmcgrccn 1999, Silk 1997, 2005) should be the topic of a future Annual 
Reviews article. 
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10 APPENDIX - NEARBY (< 1 kpc) 

MASSIVE STAR-FORMATION REGIONS 

In this Appendix, we attempt to give a complete table of star formation regions 
within Ikpc that contain at least one massive star (see Table 5). Some of 
these nearby regions of massive star formation are described in more detail by 
Bally et al. (2005). Moreover, the massive stellar content of nearby (distance 
<650pc) OB associations are disscussed in the classic Hipparcos paper by de 
Zeeuw et al. (1999), including Sco 0B2, Vel 0B2, Per 0B2, and Cep 0B2 (see 
their figures 1 and 29). We also should mention the Galactic 0-star catalog of 
Mafz-Apellaniz, Perez & Mas-Hesse (2004), which includes spectral classifica- 
tion, optical-NIR photometry, multiplicity, and astrometric information, as well 
as cluster?association membership. (See also the living 0-star Web-catalogue 



http://www-int.stsci.edu~jmaiz/GOS/GOSmain.html|. O stars in open clus- 



ter regions can also be selected from the combined photometric and astrometric 
membership analysis of 520 previously known (Kharchenko et al. 2004) and 130 
newly detected (Kharchenko et al. 2005) open clusters. As a final highlight, 
we show here Westerlund 1, the most massive young star cluster in the Galaxy 
(Figure 18). 
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Tabic 5: Star formation regions (within 1 kpc) with massive stars (B2 and ear- 
lier) 



Region 


Sources 


Distance 


Age 


L/SpT 


References 


OMC-1 


BN/KL 


450 pc 


emb" 


105 Lq 


Menten & Reid 1995 
Greenhill et al. 2004 


OMC-IS 


FIR4 


450 pc 


emb 


lO^Lo 


Schmid-Burgk ct al. 1990 
Smith et al. 2004 


NGC1977 


HD 37018 


450 pc 


5Myr 


BIV 


Makincn ct al. 1985 
Howe ct al. 1991 


NGC1980 


L Ori 


550 pc 


5Myr 


09III 


Johnstone & Bally 2006 
Piskunov et al. 2006 


NGC2023 


HD 37903 


475 pc 


~5Myr 


B1.5V 


Howe et al. 1991 
Wyrowski et al. 2000 


NGC2024 


IRS2b 


360 pc 


emb 


08V 


Bik et al. 2003 
Lenorzer et al. 2004 


cr Ori 


a Ori AB 


350 pc 


~3Myr 


09.5V 


Sanz-Forcada et al. 2004 
Hernandez et al. 2007 


Coll69 


A Ori 


400 pc 


~6Myr 


08III 


Dolan & Mathieu 2002 
Barrado y Navascues 2005 


S106 


IRS4 


600 pc 


emb 


09V 


Hodapp & Rayner 1991 

Furuya ct al. 1999 


W40 


OSla, 2a, 3a 


600 pc 


obsc'' 


09V 


Smith ct al. 1985 
Vallee & MacLeod 1994 
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Region Sources Distance Age L/SpT References 



NGC1579 


LkHalOl 


700 pc 


obsc 


B0.5V 


Barsony ct al. 1991 
Hcrbig ct al. 2004 


NGC2264 


IRSl 


800 pc 


emb 


B2V 


Schwartz et al. 1985 
Schreyer et al. 2003 


Men R2 


IRS3 


830 pc 


emb 


BIV 


Carpenter et al. 1997 
Preibisch et al. 2002 


GGD 12-15 


IRAS06084 


830 pc 


emb 


B0.5V 


Gomez et al. 1998 
Gutermuth et al. 2005 


S140 


IRSl 


900 pc 


emb 


B2V 


Lester et al. 1986 
Preibisch et al. 2001 


IC 1396 


HD 206267 


800 pc 


~3Myr 


06.5V 


Schulz et al. 1997 
de Zeeuw et al. 1999 


Cep-0B3 


Cep-A/HW2 


725 pc 


emb 


2x IO^Lq 


Garay et al. 1996 
Hartigan et al. 2000 


Cyg-X 


AFGL 2591 


1000 pc 


emb 


2x WLq 


Tamura & Yamashita 1992 
van der Tak et al. 2006 



"emb stands for embedded, i.e., a very young phase (< 1 Myr); sources optically 
invisible 

''obsc stands for obscured, i.e., a young phase (~ 1 Myr); sources optically visible 

but heavily dust cxtinctcd 

Note: AFGL 2591 may be at a distance of 1700 pc (Schneider et al. 2006) and 
part of the Cyg 0B2 association, one of the richest regions of massive stars in 
the Galaxy with ~ 100 O stars, including a hidden very massive OB cluster, 
shown in Figure 18. Note that this massive star forming region is a source 
of significant diffuse 7 ray emission, i.e. the 1.8 MeV line from the radioactive 
decay of ^^Al (Pliischke et al. 2002) 
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GLOSSARY 



• Accretion 

gas accumulation of a star or protostar, increasing the mass of the object 

• Competitive accretion 

gas initially unbound to a star or protostar and moving with relative speed, 
Vrci, w.r.t. to the object is added to it, i.e., to its gravitational sphere or 
column of influence (defined by an impact parameter 2GM/v^gi) 

• Coalescence 

growth of the mass of a stellar or protostellar object by a physical collision 
and merger with another object 



inelastic collision between two stars or protostars leading to the amalga- 
mation of the two bodies 

• Massive star 

star more massive than about 8 Mq that ends its life with a type II super- 
nova (unless it is in a close binary system with mass transfer); a massive 
star on the main sequence has spectral type B3 or earlier. 

• OB star 

star of spectral type O (> 16 Mq) or type B (B > 4 Mq) 

• Ultracompact HII region 

small blob (typical density lO^cm"^, size 0.1 pc) of ionized gas emitting 
radio continuum radiation of high emission measure (L n^ = 10^ pccm~^) 

• Hypcrcompact HII region 

very small blob (density 10^ cm^"^, typical size 0.01 pc) of ionized gas emit- 
ting radio continuum radiation of very high emission measure (L n^ — 10^'^ pc cm~^) 

• Protostar 

object on its way to become a star, with more than half of its final mass 
still to be accumulated 



Merger 
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• Infrared Dark Cloud (IRDC) 

dense interstellar cloud (often filamentary) seen in absorption against 
background thermal IR emission; detected with ISO, MSX, and Spitzer 

• Core 

basic molecular cloud unit (small dense gas fragment, mass lO-lOOM©, 
size ~0.1pc) to form one or a few stars 

• Clump 

basic molecular cloud unit (large dense gas fragment, mass ~ 1000 - 5000 Mq , 
size ~ 0.5 pc) to form a young cluster of stars 

• Hot Molecular Core (HMC) 

dense warm (T>100K) molecular gas in emission (CO, NH3, CH3OH) 
heated from a protostar inside; a compact region (size <0.1pc) in a star 
forming molecular clump such as the BN/KL region in Orion 

• Orion Nebula Cluster (ONC) 

elongated 1 x 2 pc cluster of young stars, centered on the high-mass Trapez- 
ium stars in Orion; about 4000 cluster members with a total stellar mass 
of ~1800Mq 

• Trapezium Cluster 

The 4-5 stars in the center of the Orion nebula, 9^ C, A, B, D, E in 
order of decreasing luminosity, the main exciting/ionising source being 6^ 
C (05.5V, ca. 35 M©) 

• 30 Doradus (30 Dor) 

prime example of a giant (200 pc) extragalactic HII region, powered by 
^ 100 massive O stars stars in the NGC 2070 cluster (cluster mass some 
10** M0, age <3Myr, cluster radius about 1 arcmin or 15 pc) 

• R136 

nominally the 1 arcsec core of the 30 Dor nebula and HII region (Radcliffe 
object No. 136), but R136 also often designates the central dense core of 
1 pc (4 arcsec) radius of the exciting star cluster NGC 2070 
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Jeans mass 

critical mass that must be exceeded for a gas cloud to collapse dynami- 
cally, with the self-gravity of the gas cloud being opposed only by thermal 
pressure; if squeezed by an external pressure, the critical mass for the on- 
set of collapse is called the Bonnor-Ebcrt mass. Other anisotropic forces 
such as rotation or magnetic fields can oppose the cloud's self-gravity and 
thus impede star formation. 

Free-fall timescale 

timescale it takes for a pressure- free cloud to collapse dynamically to a very 
small size (nominally a point) under its own self-gravity; the timescale de- 
pends on the average cloud gas density only (for number density n = 10^ cm" 
tfi = 105yr). 

Kelvin-Helmholtz (KH) timescale 

time it takes for a young stellar object or stellar core to radiate its thermal 
energy content, typically 10* yr for young high-mass stars but 10^ yr for 
solar-type stars (txH = GM^/RL). Nuclear energy generation eventually 
ofi^sets the thermal energy loss. 

Premain-Sequence evolution 

While high-mass stars are born on the Main-Sequence (i.e. they reach 
hydrogen burning while still accreting matter onto their stellar surface 
during their deeply embedded phase) , the collapse of a low- mass star does 
not immediately lead to a hydrogen-burning star, rather to a quasi-static 
object that keeps contracting in radius by factors of a few before hydrogen 
burning ignites in the center. During this contraction, the young stellar 
object has a time-dependent radius, effective temperature, and luminosity 
that can be calculated by pre-MS evolutionary models. 

ZAMS 

zero-age Main Sequence: a star that has contracted enough to start hy- 
drogen burning in its central region 

Initial Mass Function (IMF) 

mass distribution of (single) stars at birth introduced by Salpeter in 1955 
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• Spitzer 

NASA mid-infrared space telescope of 85 cm diameter launched in 2003 

• MSX 

mid-infrared satellite explorer with an aperture of 20 cm diameter 

• Chandra 

0.5-10 kcV NASA X-ray telescope in space with particularly good imaging 
quality comparable to optical seeing (angular resolution 1 arcsec) 

• HST 

Hubble Space Telescope, 2.4 m aperture, workhorse for high angular reso- 
lution space astronomy in UV, optical, and near-IR; in orbit since 1990 
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Figure 1: Accretion and mass loss as exchange between components: the ac- 
cretion disk as reservoir and interface between the molecular cloud core and the 
forming star. 
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Figure 2: Hubble Space Telescope optical/IR image of the dense mas- 
sive young cluster R136/30 Dor (courtesy of M.J. McCaughrean; FOV 
~30arcsec x 30arcsec or 7.5 pc x 7.5 pc). Dozens of massive O stars are found 
crowded within the half-light radius of 2pc (Brandl et al. 1996). (a) A VLT 
image of NGC 3603 (Brandl et al. 1999) and (b) a VLT image of the Trapezium 
Cluster in Orion (McCaughrean 2001) are shown, as these two galactic clusters 
would be seen if they were located at the distance of R136 in the Large Magel- 
lanic Cloud (50kpc) and imaged with similar angular resolution (see Zinnecker 
2002). 
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Figure 3: Collapse and fragmentation of a giant molecular cloud, simulated 
with smoothed particle hydrodynamics (SPH) using sink particles (LA. Bon- 
nell, P.C. Clark & H. Zinnecker, in preparation). For a description of SPH, see 
the box. Plotted is the spatial distribution of gas column density, color-coded 
such that deep blue refers to the highest values. The initial conditions for this 
three-dimensional simulation included a molecular gas cloud of mass 10^ Mq and 
diameter 100 pc, somewhat centrally condensed (factor of 20, with a Gaussian 
radial profile). The cloud's turbulent kinetic energy was equal to its gravita- 
tional energy. Note the filamentary structure and the associated dense cores 
(blue dots). These cores technically are represented by sink particles (~2400 
at this stage of evolution). With typical masses of 10 to 100 M0 and typical 
sizes of 0.1 pc they could be the initial fragments for high- mass star formation. 
A total of 2.5 million SPH particles was used. Some similarity to cosmological 
simulations of structure formation is noted. 
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Figure 4: ZAMS luminosity (blue), ZAMS radius (red dashed), and Kelvin- 
Helmholtz quasi-static gravitational contraction timescale toward the ZAMS 
(green dotted) as a function of stellar mass (ZAMS means zero-age main se- 
quence) . These values are extracted from t = yr models given by Meynet & 
Maeder (2005) for rotating stars M*/M0 > 12 with mass loss, by Pietrinferni et 
al. (2004) for stars 0.5 <M*/M0 < 10, and for completeness from unpublished 
tracks by Yorke for masses M» =0.1, 0.2, and 0.3 Mq, using the computer code 
provided in the book by Bodenheimer et al. (2007) . ZAMS radii are well repre- 
sented by two power laws with a break at M^/M© = 1.5. The power law slope 
of the luminosity-mass relation for massive stars varies from about 3.7 to 1.6 in 
the range 8 < MJMq < 120. 
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Figure 5: ZAMS luminosity to mass ratio (left scale) and corresponding critical 
effective opacity (right scale) as defined by Equation 3 as a function of stellar 
mass. The solid grey line depicts the ZAMS models shown in Figure 4. Dotted 
grey lines indicate the total luminosity (including accretion luminosity) of stars 
accreting at the indicated constant rate as discussed in the text (cf. Figure 8). 
The triangle denotes the position of an 05V star (see Figure 9). Dashed 
lines denote the opacities of dusty gas in the light of black body sources at the 
temperature indicated (for two different grain types, see Figure 6) and the 
contribution from electron scattering in a fully ionized plasma. Because of the 
combined effect of UV lines (UV line blanketing), the actual opacity of a hot 
plasma can be greater than that from electron scattering alone. This may be 
relevant for the existence of an upper stellar mass limit. 
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Figure 6: Planck- weighted mean effective opacity of dusty gas, using the 
(Preibisch et al. 1993) dust model with ice-coated grains (blue curve) and 
grains without ice mantels (red curve), assuming black-body radiation at the 
temperature Tj-ad and solar abundances. The right-hand scale (grey) is based 
on the ZAMS models shown in Figure 5 (solid grey line). 
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Figure 7: Edge-on cut through an accreting protostar and its circumstellar disk 
and envelope. The protostar is located at (0,0) and is too tiny to see. Tem- 
perature (red scale, upper-left), density (blue scale, upper-right), and velocity 
(length and direction of arrows) distributions of the accreting material are dis- 
played (case F of Yorke & Bodenheimer 1999, recalculated using an improved 
version of their code). White contour lines are plotted for log T = 2.0, 2.5, and 
3.0 and for log/? = -18, -16, -14, and -12. At this evolutionary time t — 65,000 yr 
after formation of the protostellar core, 7.0 M0 of material have accreted onto 
the protostar, 2.8 are in the disk, and 0.2 are in the infalling envelope. 
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Figure 8: Evolutionary tracks in the HR diagram for (proto-)stars accreting at 
a constant rate (colored lines) are contrasted to the tracks of non-accreting stars 
(black lines). All accreting tracks are assumed to begin at the birthline of an 
equilibrium deuterium burning 0.1 pre-main-sequence star. Non-accreting 
tracks up to H burning were calculated by Yorke for M^/M© = 0.1, 0.3, 1, 3, 10, 
and 30, using the computer code supplied with the book by Bodenheimer et al. 
(2007). [Adapted from Yorke (2002)] 




Figure 9: Luminosity and effective temperature calibration of luminosity classes 
f (blue- filled circles), fff (green squares), and V (red triangles) for stars of 
spectral classes 03, 04, 05, 05.5, 06, 06.5, 07, 07.5, 08, 08.5, 09, and 09.5. 
For comparison, theoretical evolutionary tracks from Meynet & Maeder (2003) 
for non-rotating (dashed grey) and for rotating stars (continuous grey) with mass 
loss are labeled by their starting ZAMS mass. Meynet et al. (1994) isochrones 
for non-rotating stars are also plotted for (solid black line labeled ZAMS), 1, 
3, and 5 Myr (dotted black lines). [Adapted from Figure 14 of Martins, Schaerer 
& HiUier (2005)] 
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Figure 10: Flux of radiation that can ionize hydrogen (blue lines) and neutral 
helium (red dashed lines) as a function of age (counted from the moment of ar- 
rival on the ZAMS), for stars of ZAMS masses (top to bottom) of M^/M© = 120, 
85, 60, 40, 25, and 20. Mass loss is included in the evolution. [Based on models 
from Schaerer & de Koter (1997)] 
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Figure 11: The inner accretion disk around a close massive accreting binary 
pair: inward radial flow is allowed in the equatorial plane. A polar cavity is 
evacuated by a combination of radiation and the stellar wind. The disk is 
self-shielded from the intense EUV fleld by an ionization front separating HII 
and HI gas; and from the FUV fleld interior to the HI/H2 interface by dust, 
molecular hydrogen, and CO. The dust is destroyed at rdust- Interior to rcvap, 
the radius where the sound speed exceeds the escape velocity, even the ionized 
gas is gravitationally bound. Sizes are not to scale. 
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Figure 12: Time-dependent evolution of a turbulent, self-gravitating 1000-Mq 
cloud showing the formation of sheets, filaments, and cores, the latter of which 
become gravitationally unstable. Newly formed stars are shown in dark blue, 
the gas is shown in blue-green. Stars tend to cluster, they continue to accrete 
material in competition with other stars. Times shown are 1.0, 1.4, 1.8, and 2.4 
initial free-fall times (tff), from left to right and top to bottom. [Adapted from 
Bonnell, Bate & Vine (2003)]. 
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Figure 13: (a) A grazing encounter between two 3M0 premain-sequence stars 
that results in the formation of a binary (ri„in = 25.8 solar radii). (&) A de- 
tached encounter between a 3M0 premain-sequence star and a 10 M0 zero age 
main-sequence star with the same minimum periastron distance as in the (a) 
encounter. The stars have radii of 12.9 and 3.92 solar radii, respectively. The 
greater density of the 10 Mq star results in tidal disruption of the low-density 
3 Mq star to form a disk around the massive star. The encounters have zero 
relative velocity at infinity (i.e., they are parabolic encounters). [ Adapted from 
Zinnecker & Bate (2002), Davies et al. (2006)]. The cross section for a subse- 
quent collision with (say) another 10 Mq star is significantly increased, and a 
runaway coUisional growth in mass is possible. This is the so-called 'shred and 
add' process. 
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Figure 14: Massive stars (dark blue circles) are formed in the center of individ- 
ual subclusters of low-mass stars (light circles) because of competitive accretion. 
These subclusters evolve by merging with the other subclusters. The final state 
of the simulation is a single, centrally condensed cluster with little substructure 
but with 4 massive stars, one from each subcluster. This then is a model for the 
origin of Trapezium-type systems. [Adapted from Bonnell, Vine & Bate (2004)] 
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Figure 15: The multiple star systems of the Orion Trapezium as revealed by 
bispectrum speckle interferometry. [Courtesy of G. Weigelt & Th. Preibisch; 
inserts adapted from Schertl et al. (2003) and Kraus et al. (2007)] 
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Figure 16: A wide-field JHKs composite image of the AFGL 2591 massive 
outflow taken with the NIRI camera at the Gemini North telescope in excellent 
seeing (courtesy of C. Aspin and Gemini Observatory; FOV ~2arcmin, see- 
ing 0.35 arcsec). The broad-band Ks filter includes the ro- vibrational molecular 
hydrogen v = 1-0 S(l) emission line at 2.12 /xm, which is indicative of shock- 
excited gas. However, much of the nebulosity is likely seen in reflection because 
it is also detected in J and H. Note the multiple poorly-coUimated flows, loops, 
and cavities emanating to the west from this massive (ca. 20 M©) star (see also 
the speckle image in Preibisch et al. 2003 from the 6m SAO telescope). Nothing 
of this kind is seen to the east of the source, except for a few faint bow-shocks, as 
the (redshifted) counterflow is deeply embedded and probably partly obscured 
by a dense circumstellar disk (Trinidad et al. 2003) around this very young 
star. The fact that the flow and counterflow both reach the same distance of 
'--^ 40 arcsec (0.2 pc) from the central source indicates that this is intrinsically a 
bipolar outflow, confirmed by high-velocity CO radio observations (Lada et al. 
1984; see also the H2 and optical spectroscopic observations of the associated 
Herbig-Haro objects by Tamura & Yamashita 1992 and Poetzel, Mundt & Ray 
1992, respectively). It seems that the bipolar outflow from this late 0-type 
massive star has both a well-coUimated and a wide-angle component! 
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Figure 17: Schematic representation of a self-similar fractal hierarchy of dense 
gas cores in a turbulent molecular cloud with three branches per level (courtesy 
J. Melnick). The chance of randomly picking a core of mass M in a cloud with 
such hierarchical structure is proportional to 1 /M. Because the cloud branches 
correspond to log intervals, random sampling yields a mass spectrum of (prestel- 
lar) cores with equal mass in equal logarithmic intervals, i.e., a power-law 
dN/dlogM with index -1. If there is competition for mass, i.e., if some smaller 
(and presumably denser) cores are turned into stars before the bigger cores, of 
which they are part, have time to collapse, then there will no longer be enough 
mass available to form a star with the mass of the undiminished larger core. 
Thus, high-mass cores (stars) get depleted from the mass function if low-mass 
cores (stars) form first. The faster formation of low- mass stars, combined with 
the mass depletion of subparts inside clouds, converts the power-law mass spec- 
trum with a power index of -1 into one with an index of -1.35 (or so), i.e., into a 
Salpeter initial mass function. Note that in this reasoning (following Elmegreen 
1997) it was assumed that a fixed fraction of each and every core mass is turned 
into stellar mass (single, binary, or multiple). 
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Figure 18: Composite colour JHKs image of the ~4Myr old Galactic starburst 
cluster Westerlund 1 (Wd 1) in Cyg 0B2 obtained with the ESQ NTT telescope 
(courtesy W. Brandner). This obscured object is held to be the most massive 
young cluster in the Galaxy, with ca. 100 0-type stars and a mass of at least 
20,000 Mq (see Knodelseder 2000, Clark et al. 2005, Mengel & Tacconi-Garman 
2007, Brandner et al. 2007). A number of massive eclipsing binary systems have 
been identified (Bonanos 2007) as well as a very massive runaway star, an Of- 
type supergiant, whose origin can be traced back to the cluster (Comeron & 
PasquaU 2007). Like the Orion Nebula Cluster, or the NGC 3603 and 30 Dor 
central clusters, Wd 1 is strongly mass segregated with the high mass stars 
being more centrally concentrated than the low mass stars. The field of view is 
4' X 4', corresponding to 4.2 pcx 4.2 pc at the distance of Wd 1. 



